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ABSTRACT 

We study the physical conditions in damped Lyman-a systems (DLAs) , using a sample 
of 33 systems toward 26 QSOs acquired for a recently completed survey of H2 by 
Ledoux et al. (2003). We use the column densities of H2 in different rotational levels, 
together with those of C i, C i*, C i**, C 11* and singly ionized atomic species to 
discuss the kinetic temperature, the density of hydrogen and the electronic density 
in the gas together with the ambient UV radiation field. Detailed comparisons are 
made between the observed properties in DLAs, the interstellar medium (ISM) of the 
Galaxy, the large and small Magellanic clouds (LMC and SMC). 

The mean kinetic temperature of the gas corresponding to DLA subcomponents 
in which H2 absorption line is detected, derived from the ortho-to-para ratio (153±78 
K), is higher than that measured in the ISM (77±17 K) and the Magellanic clouds 
(82±21 K). Typical pressure in these components (corresponding to T = 100—300 
K and tt-h = 10—200 cm~^), measured using C i fine-structure excitation, are higher 
than what is measured along ISM sightlines. This is consistent with the corresponding 
higher values for A''(H2,J=2)/iV(H2,J=0) seen in DLAs. From the column densities of 
the high-J rotational levels, we derive that the typical radiation field in the II2 bearing 
components is of the order of or slightly higher than the mean UV field in the Galactic 
ISM. Determination of electron density in the gas with II2 and C i show the ionization 
rate is similar to that of a cold neutral medium (CNM) in a moderate radiation field. 
This, together with the fact that we see II2 in 13-20% of the DLAs, can be used to 
conclude that DLAs at z > 1.9 could contribute as much as 50% star formation rate 
density seen in Lyman break galaxies (LBGs). 

C II* absorption line is detected in all the components where H2 absorption line is 
seen. The excitation of C il in these systems is consistent with the physical parameters 
derived from the excitation of II2 and C i. We detect C 11* in about 50% of the DLAs 
and therefore in a considerable fraction of DLAs that do not show II2. In part of 
the later systems, physical conditions could be similar to that in the CNM gas of 
the Galaxy. However, the absence of C i absorption line and the presence of Al ill 
absorption lines with a profile similar to the profiles of singly ionized species suggest 
an appreciable contribution from warm (WNM) and/or partially ionized gas. The 
absence of H2, for the level of metallicity and dust depletion seen in these systems, 
are consistent with low densities (i.e njj ^ 1 cm~'^) for a radiation field similar to the 
mean Galactic UV field. 

Key words: Cosmology: observations - Galaxies: halos - Galaxies: ISM - Quasars: 
absorption lines 
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1 INTRODUCTION 

Damped Ly-a (DLA) systems seen in QSO spectra are char- 
acterized by very large neutral hydrogen column densities: 
7V(H i) > 2 X 10^" cm~^. Such an amount of neutral gas 
is usually measured through local spiral disks. The case for 
DLA systems to arise through proto-galactic disks is further 
supported by the fact that the cosmological density of the 
absorbing gas at Zahs ~ 3 is of the same order of magni- 
tude as the cosmological density of stars at present epochs 
(Wolfe 1995). Moreover, the presence of heavy elements 
{Z ~ 1/10^0 ) suggests that DLAs are located in over-dense 
regions where star formation activity takes place (Pettini et 
al. 1997) and at low and intermediate redshifts strong metal 
line systems and DLAs have been demonstrated to be asso- 
ciated with galaxies (e.g. Bergeron & Boisse, 1991; Le Brun 
et al. 1997). It has also been shown that the profiles of the 
lines arising in the neutral gas show evidence for rotation 
(e.g. Prochaska & Wolfe 1997). However, hydrodynamical 
simulations have shown that the high rodshift progenitors 
of present-day galactic disks could look like an aggregate of 
well separated dense clumps. In fact, the kinematics seen in 
the absorption line profiles of DLAs could be explained by 
relative motions of the clumps with little rotation (Haehnelt 
et al. 1998; Ledoux et al. 1998). 

Studying the star-formation activities in DLAs is very im- 
portant for the understanding of galaxy formation in the 
Universe. Recently, Wolfe et al. (2003a, 2003b, 2004) have 
shown that, even if DLAs sustain only a moderate star- 
formation activity, they will contribute appreciably to the 
global star-formation rate (SFR) density at high redshifts. 
The SFR in DLAs can be estimated either by detecting the 
galaxies responsible for DLAs or by inferring the intensity of 
the UV field in DLAs using the induced excitation of atomic 
and molecular species. In the latter case, it is important to 
have a clear understanding of the physical conditions in the 
gas to derive an accurate estimate of the SFR. In the case of 
the Galactic ISM, rotational excitations of H2 (see Brown- 
ing et al. 2002 and references there in) and fine-structure 
excitations of C I, C 11, O I and Si 11 are used to derive the 
physical state of the absorbing gas (see for example Wclty et 
al. 1999). Detecting and studying these transitions in DLAs 
is the first step toward understanding the physical condi- 
tions and hence the star-formation activity in DLAs. 
Molecular hydrogen is ubiquitous in the neutral phase of the 
interstellar medium (ISM) of galaxies. Formation of H2 is ex- 
pected on the surface of dust grains, if the gas is cool, dense 
and mostly neutral, and from the formation of H~ ions if the 
gas is warm and dust-froo (see e.g. Jenkins & Peimbert 1997; 
Cazaux & Tiolens 2002). As the former process is most likely 
dominant in the neutral gas associated with DLA systems, 
it is possible to obtain an indirect indication of the dust 
content in DLAs without depending on extinction and/or 
heavy element depletion effects. Moreover, by determining 
the populations of difi'erent H2 rotational levels, it is possi- 
ble to constrain kinetic and rotational excitation tempera- 
tures and particle densities. Effective photo-dissociation of 
H2 takes place in the energy range 11.1 — 13.6 eV through 
Lyman- and Werner-band absorption lines and the intensity 
of the local UV radiation field can therefore be derived from 
the observed molecular fraction. A direct determination of 
the local radiation field could have important implications in 



bridging the link between DLA systems and star-formation 
activity at high redshifts. 

We have searched for molecular hydrogen in DLA and sub- 
DLA systems at high redshift {zahB > 1.8), using UVES 
at the VLT down to a detection limit of typically A'^(H2) ~ 
2x 10" cm-2 (see Ledoux et al. 2003). Out of the 33 systems 
in our sample, 8 have firm and 2 have tentative detections 
of associated H2 absorption lines. In all of the systems, we 
measured metallicities relative to Solar, [X/H] (with either 
X=Zn, or S, or Si), and depletion factors of iron, [X/Fo], 
supposedly onto dust grains. Although H2 molecules are de- 
tected in systems with depletion factor, [Zn/Fe], as low as 
0.3, the systems where H2 is detected are usually amongst 
those with the highest metallicities and depletion factors. 
In particular, H2 is detected in the three systems with the 
largest depletion factors. Moreover, in two different systems, 
one of the H2-detected components has [Zn/Fe] > 1.5. This 
directly demonstrates that a large amount of dust is present 
in the components where H2 is detected. The moan H2 
molecular fraction, / = 2iV(H2)/[2iV(H2) -hiV(Hl)], in DLA 
systems is generally small (typically log / < —1) and sim- 
ilar to what is observed in the Magellanic Clouds. There 
is no correlation between the amount of molecules and the 
neutral hydrogen column density; in particular, two systems 
where H2 is detected have log A'^(Hl) < 20.3. Approximately 
50 percent of the systems have log / < —6: this is probably 
a consequence of a reduced formation rate of H2 onto dust 
grains (probably because the gas is warm, T > 1000 K) 
and/or of an enhanced ionizing flux relative to what is ob- 
served in our Galaxy. 

In this work, we present additional high S/N ratio data on 
three of the DLA systems in which H2 is detected and the 
results of multi-component Voigt profile fits to neutral and 
singly ionized species (including C I, C I* and C 11*) in all 
the DLAs in our sample. We estimate the range of physi- 
cal conditions in the neutral gas using standard techniques 
that are used in ISM studies. The paper is organized as fol- 
lows. In Section 2, wo give the details of the additional data 
and present the new fits to the H2 absorption lines in the 
corresponding three systems. In Section 3, we discuss the 
relative populations of different H2 rotational levels deriv- 
ing information on the physical state of the gas by com- 
paring the DLA observations with Galactic ISM, SMC, and 
LMC data. In Sections 4 and 5 we discuss, respectively, the 
fine-structure excitation of C I and the ionization state of 
Carbon. In Section 6 we study the C 11* excitation in de- 
tail. Finally, we summarize our results and discuss various 
implications of the overall study in Section 7. 



2 DATA SAMPLE 

The Ultraviolet and Visible Echelle Spectrograph (UVES; 
Dokkor et al. 2000), installed at the ESO VLT 8.2-m tele- 
scope unit Kueyen on Mount Paranal in Chile was used to 
search for H2 in a large sample of DLAs. The sample and 
data reduction procedure are described in detail in Ledoux 
et al. (2003). Observations and details of Voigt profile analy- 
sis of H2 and metal line absorption linos toward Q 0013—004, 
Q 0551-366 and Q 1232-1-082, along the lines of sight of 
which H2 is detected, are described in, respectively, Petit- 
jean et al. (2002), Ledoux et al. (2002) and Srianand et al. 



UVES-VLT survey for H2 in DLA systems 3 



L 14-0 & W 3-0 


/ ^ If / f V pf f K 

V^- K,,;iiv-\p\.-^^^ 1 


ft W i 


3405 

L 9-0 


3410 34 15 3420 


1 I I i 








3565 

L 5-0 

1 1 1 1 1 1 


3570 3575 3560 


3565 


- -. . 1 . 


Vy~-U....™p/: - - _ ,- W-j™-..-:^- - - , , 1 , , , p 


■ r 


3730 

L 3-0 


3735 3740 3745 


3750 




4000 

Observed wavelength (A) 



W 2-0 & L 12-0 & L 11-0 




Figure 1. Voigt profile fits to different rotational levels in the ground vibrational state of H2 detected at Zabs= 2.5948 toward Q 0405—443 
(top panels) and 2abs= 2.8111 toward Q 0528— 250 (bottom panels). The fit parameters are given in Table [Tl 



(2000). The Voigt profile fits to H2 and other metal lines 
at Zabs= 3.024 toward Q0347 - 383, Zabs= 2.595 toward 
Q 0405 - 443 and Zabs= 2.0868 toward Q 1444 + 014 are dis- 
cussed in Ledoux et al. (2003). For systems in which H2 is 
not detected, Ledoux et al. (2003) have provided upper lim- 



its on Af(H2) together with mean metallicities and depletion 
factors. 

Recently, we have obtained additional higher spec- 
tral resolution spectra {R ~ 55,000) of Q 0347—383 and 
Q 0405 — 443 as a part of our ongoing programme on cos- 
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Table 1. Results of Voigt profile fitting to H2 absorption lines 
using the new data 



QSO 



2abs J log7V{cm 2) fe (km s ^) 



Q 0405 -443 2.59471 



2.59486 



Q 0528 - 250 2.81100 



2.81112 



Q 1232 + 082 2.33772 



14.81^ 



17.65lH^ 

17.97in^8 
15 93+0 ''5 
1+0.14 
-0.26 

< 13.55 

15.24lg;J3 

14.00ini 
10 Qo+0.09 

10.30_g U5 

< 13.55 
,7 29+0 08 
^'■^^-0.16 

17.78_o 

15.68lg:f, 

14.13in} 
13.681H1 
17.24lHf 
17.68ini 
16.98tg;J^ 

in c:o + 0.02 

lCi.Oi5_Q Q4 

10 34+0.10 

19 17+010 
^"■^'-0.10 

16.70t°i2 
16.90t°j2 

14.68in^5 
14.451H1 



1.5±0.5 



l.litO.l 



2.8±0.5 



1.0±0.3 



4.5±0.5 



mic variation of the electron-to-proton mass ratio (Petitjean 
et al. 2004). Nine exposures of 1.5 h each were taken for 
each of the quasars over six nights under sub-arcsec seeing 
conditions in January 2002 and 2003 for Q 0347—383 and 
Q 0405 — 443, respectively. We have also obtained additional 
data of Q 1232-1-082 to study the HD lines that are detected 
in the DLA (Varshalovich et al. 2002). Spectra were reduced 
using the UVES pipeline and addition of individual expo- 
sures were performed using a sliding window and weight- 
ing the signal by the errors in each pixel. We detect a new 
H2 component at Zabs~ 2.59486 toward Q 0405—443 in ad- 
dition to the strong component reported in Ledoux et al. 
(2003). We also present Voigt profile fits to the H2 lines in 
the Zabs= 2.811 system toward Q0528 - 250. The single H2 
component seen in the lower spectral resolution CASPEC 
spectrum (Srianand & Petitjean 1998) is resolved into two 
distinct components in our new UVES spectra. For both 
these systems the Voigt profile fits to the H2 Lyman and 
Werner band absorption lines are shown in Fig. Q and re- 
sulting parameters are summarised in Table This Table 
also gives the results of Voigt profile fits to H2 for the 2abs= 
2.33772 toward Q 1232+082. 

The main purpose of this paper is to provide a de- 
tail account of C I, C 11* and other metal lines in DLAs 
of our sample and extract physical conditions in conjunc- 
tion with the H2 content reported in Ledoux et al. (2003). 
For the 2abs= 2.139 system toward Tol 10374-014 and the 



-Zabs= 3.350 system toward Q1117 — 1329, we use the re- 
sults presented in Srianand & Petitjean (2001) and Peroux 
et al. (2002) respectively. For the rest of the systems we 
give here the results of the multicomponent Voigt profile 
fits. For this we use a Voigt-profile fitting code that de- 
termines the best fitting parameters (column density, ve- 
locity dispersion and redshift) using minimization tech- 
niques (Chand et al. 2004). We use the oscillator strengths 
compiled in Table 1 of Ledoux et al. (2003) for metal ions 
and those given by Morton & Dinerstein (1976) for H2. 
In this article, we measure metallicities relative to Solar, 
[X/H] = log[iV(X)/iV(H)] - log[iV(X)/iV(H)]Q, with either 
X=Zn, or S, or Si, and depletion factors of iron, [X/Fe] = 
log[iV(X)/iV(Fe)] - log[iV(X)/7V(Fe)]0, adopting the Solar 
abundances from Savage & Sembach (1996). 



3 DETERMINATION OF PHYSICAL 
PARAMETERS USING H2 LEVEL 
POPULATION 

In this section, we estimate different physical parameters 
from the column densities of H2 in different J rotational 
levels. 



3.1 Kinetic temperature of the gas 

It is a standard procedure, in ISM studies, to use the ortho- 
to-para ratio (OPR) to infer the kinetic temperature of the 
gas assuming local thermodynamic equilibrium, LTE (Tum- 
linson et al. 2002, and references there in). Indeed, recent 
numerical investigations suggest that the OPR is a good 
tracer of the kinetic temperature over large regions of the 
parameter space (Shaw et al. 2004). For completeness, we 
first review our understanding of the OPR and outline the 
method for deriving the kinetic temperature before applying 
the method to the data. 



3.1.1 General outline 

As the interconversion between para and ortho states in- 
volves a spin flip, it is not allowed for processes involving an 
isolated molecule (i.e., radiative processes cannot induce in- 
terconversion). Ortho/para interconversion is only possible 
through (i) spin exchange induced by collisions with protons 
(with a rate coefficient in the range 10~^° - 10"'' cm s ; 
see Dalgarno, Black & Weisheit 1973; Flower & Watt 1984 
and Gerlich 1990) or with hydrogen atoms (with a rate co- 
efficient an order of magnitude less than that of protons; 
Mandy & Martin 1993; Tine et al. 1997) and (ii) reactions 
on the surface of dust grains (Le Bourlot 2000). In the case 
of local thermodynamic equilibrium (LTE) , 



OPRlte = 3 



^j^„,,(2J + l) exp [-BJ(J + 1)/T1 
Ej=ovcn(2J + l) exp [-BJ(J + 1)/T] 



(1) 



where, J, is the rotational quantum number, B is the rota- 
tional constant of II2 {B — 85.3 K), and T is either the ki- 
netic temperature of the gas (when OPR is governed by spin- 
exchange collisions) or the formation temperature (when 
OPR is governed by H2 formation on the surface of dust 
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Table 2. Excitation temperatures measured in individual H2 components of DLAs in the sample of Ledoux et al. (2003) 
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7gg+584 
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16 
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20 
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18 


1 c: + 0.15 
^"•-0.15 
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Note: The 2^1,5= 1.97296 and 1.97380 components toward Q 0013— 004 are not considered here as the H2 column densities 
of the J = and J = 1 levels are not known accurately due to saturation effects (see Petitjean et al. 2002) 
" Integrated H I column density in all the components (with and without H2). 
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Figure 2. The Ortho-to-Para ratio (OPR) obtained using the 
populations of the 5 rotational levels in individual components 
of DLAs (black dots) is plotted against the total H2 column den- 
sity. Other data points are from Savage et al. (1977), Spitzer, 
Cochran & Hirshfelf (1974) for the Galactic ISM (triangles), and 
Tumlinson et al. (2002) for the LMC (squares) and SMC (aster- 
isks). The vertical short-dashed line shows the high temperature 
LTE limit of the OPR (i.e., OPR = 3). 



Figure 3. The rotational temperature Toi obtained from the 
A'^(J=l)/Af(J=0) ratio in DLAs (black dots) is plotted against the 
total column density of molecular hydrogen. Other data points are 
from Savage et al. (1977), Spitzer, Cochran & Hirshfelf (1974) for 
the Galactic ISM (triangles), and Tumlinson et al. (2002) for the 
LMC (squares) and SMC (asterisks). The vertical short-dashed 
lines show the mean and Icr range of Tgi measured by Savage et 
al. (1977) in the Galactic ISM. 



grains -with LTE distribution characterized by the forma- 
tion temperature Tform; see Sternberg & Neufeld 1999; Taka- 
hashi 2001). The equilibrium temperature, T(OPR), can be 
obtained using the observed value of OPR and Eq. This 



will trace the kinetic temperature of the gas if spin-exchange 
collisions are mainly responsible for the observed OPR. 

If the gas is dense and cold and if most of the 
H2 molecules are in the J = and J = 1 levels then, 
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OPRlte 



N(J ^ 1) 
N(J = 0) 



= 9 X exp(-170.5/To 



(2) 



In the case of very optically thick molecular gas for which 
there is enough self-shielding, the A'^(J=1)/A''(J=0) ratio can 
be maintained at its Boltzmann value and the excitation 
temperature, Tqi, equals the kinetic temperature. Savage et 
al. (1977) measured a mean excitation temperature, Toi = 
77±17 K, for the galactic ISM. This is consistent with the 
mean temperature of the ISM measured using 21 cm ab- 
sorption lines. Thus it is widely believed that when there is 
sufficient shielding (i.e log A''(H2) cm~^ J5 16.5), Toi is a rea- 
sonably good tracer of the kinetic temperature. This is be- 
cause in the shielded region, II2 photodissociation time-scale 
can be larger than the time-scale for charge exchange colli- 
sion (Flower & Watt, 1984). Also a recent multi- wavelength 
study of Galactic sightlines show the Toi measured in op- 
tically think cases closely follow the spin temperature mea- 
sured from 21 cm observations (see Roy et al 2005). 

The excitation temperature, Tij , between different rota- 
tional levels (say J= i and j) of a given species (either ortho 
or para H2) can be obtained using, 



N{J ■ 



2j + l 



exp(-B[j(j + l)-*(* + l)]/T..). (3) 



N{J^i) 2i + l 

Unlike OPR, this ratio can be altered by radiation pumping 
and formation pumping in addition to collisions. If collisions 
dominate the rotational excitation then Tij will be equal to 
T(OPR). Presence of formation pumping and/or UV pump- 
ing will make Tij > T{OPR). In the following section we 
discuss various temperature estimates from the DLA sam- 
ple. 



3.1.2 Kinetic temperature of the H2 components 

In our sample, H2 is only detected in JsJ 5 levels of the vibra- 
tional ground state. Thus we compute the OPR by summing 
the II2 column densities for levels with J ^ 5. The observed 
value of the OPR for each DLA is given in column #8 of 
TableU We calculate T(OPR) from the measured OPR for 
individual systems using Eq. 0^see column #9 of Table 
When the kinetic temperature (or formation temperature) is 
high (i.e., T ^ 200 K) the OPR reaches 3, the value expected 
based on spin statistics. For a kinetic temperature similar 
to that seen in the cold neutral medium of our Galaxy (~ 
80 K) the expected OPR under LTE assumption is less than 
1. From Table |5| it is clear that the LTE temperatures mea- 
sured from the OPR for DLAs are most of the times higher 
than 80 K (the mean found in the Galactic ISM.) 

In Fig.|5|we plot the observed values of the OPR against 
the total H2 column density in the ISM (triangles), LMC 
(squares), SMC (asterisks) and DLAs (circles with error- 
bars). It is apparent that most of the OPR values in DLAs 
are significantly different from 3 (see also column #8 of Ta- 
bleHJ. The distribution of the OPR as a function of iV(Il2) 
in DLAs is consistent with that observed along Galaxy, LMC 
and SMC sightlines (see Fig.|2l when log N(H2)<16 and for 
rest of the components OPR in DLAs are systematically 
higher than that measured in the Galaxy, LMC and SMC. 
For example, OPR ^ 3 is seen only along sightlines with low 
II2 optical depth (i.e A''(H2) ^ 10^** cm"^) in the Galaxy, 
LMC and SMC (see Fig.|lJ. On the contrary, out of the three 



DLA components with OPR ^ 3, two, at 2:abs= 2.08680 and 
Zshs— 2.08692 toward Q 1444-1-014, are optically thick in the 
Lyman band absorption lines. 

We next investigate the dependence of Toi (measured 
using Eq. |5J on the total H2 column density (see Fig|^. 
Individual values measured in DLAs are listed in Table |21 
(see column #5). The large errors on both A'"(H2) and Toi 
are mostly a consequence of the difficulty to measure the 
Doppler parameter when the lines are saturated. In the case 
of the Zabs= 1.96685 component toward Q0013 — 004 the 
uncertainty is a consequence of line blending (see Petitjean 
et al. 2002). The vertical dotted lines show the mean and 
1(7 range of Toi measured by Savage et al. (1977). The data 
points from the Magellanic clouds (Tumlinson et al. 2002) 
are consistent with this range (mean Toi = 82±21 K). As in 
the case of the OPR, most of the measurements from DLAs 
with optically thick H2 (i.e log N(H2)^ 16.5) are well sepa- 
rated from that of the ISM and Magellanic clouds (Fig. 
and the spread seen in the optically thin case is consistent 
with that seen in local ISM. Note that the system with 
lowest molecular content (2abs= 3.02489 toward 0347 — 383) 
has A^(J=0) an order of magnitude lower than A'^(J =1). Toi 
can not be computed in this case as the maximum expected 
column density ratio, A'^(J=1)/A'^(J=2), is 9 under LTE con- 
ditions (see Eq. |5J. For the high optical depth clouds (i.e, 
log A''(H2) cm"2 > 16.5) in DLAs the mean Tqi is 153±78 
K. In most of the components the two temperatures Toi and 
T(OPR) are consistent within errors. This is mainly because 
most of the H2 molecules reside in the ground states. 

In summary, if we assume LTE then T(OPR) and Toi 
measured in DLAs (with log N(H2)S5l6.5) at high redshift 
are on an average higher than that measured in the ISM, 
LMC and SMC sightlines. In this high N(H2) range Tqi is 
expected to trace the kinetic temperature. However, in the 
case of optically thin systems the T(OPR) (or Toi) measured 
in DLAs are consistent with that measured in LMC, SMC 
and Galactic sightlines. Under the LTE assumption we find 
that II2 components in DLAs have kinetic temperatures in 
the range 100-200 K. 

3.2 Rotational excitation 

The rotational level populations are affected by particle col- 
lisions, UV pumping, and formation pumping. While the 
coUisional excitation plays a significant role in populating 
the low-J levels, those with Jj53 are usually populated by 
formation processes and UV pumping. In what follows, we 
discuss the excitation of H2 as seen in DLAs and compare 
with ISM, LMC and SMC sightlines. 

3.2.1 Low-J excitation 

The coUisional contribution to the excitation of H2 can 
be investigated by studying the A'^(J=2)/A'^(J=0) and 
Af(J=3)/Af(J=l) ratios. In general J = 2 and J = 3 lev- 
els can also be populated by deexcitation of H2 formed in 
the high-J states (usually referred to as formation pumping) 
or through UV pumping. The coUisional excitation rate for 
the J=0^2 transition is about an order of magnitude higher 
than that of the J=l— >3 transition for kinetic temperatures 
in the range 100 to 300 K (Forrey et al. 1997). The sponta- 
neous decay rate from J = 3 is an order of magnitude smaller 
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Figure 4. Critical hydrogen density {njj) as a function of kinetic 
temperature for the thermalization of the J = and 2 levels. The 
curve is obtained by equating the coUisional de-excitation and 
spontaneous decay rates. We have used the coUisional rates at 
low temperature given by Forrey et al. (1997). The spontaneous 
decay rates are from Turner et al. (1977). It is clear from the 
short-dashed lines that for lOOsCTsCSOO K the critical density is 
175^nH (cm-3)^60. 



Figure 5. The A''(J=2)/A''(J=0) ratio observed in individual 
DLA components (black dots) is plotted against the total molec- 
ular hydrogen column density. Other data points are from Savage 
et al. (1977), Spitzer, Cochran & Hirshfelf (1974) for the Galac- 
tic ISM (triangles), and Tumhnson et al. (2002) for the LMC 
(squares) and SMC (asterisks). The vertical dotted lines are the 
expected values of the column density ratio assuming the LTE for 
four different excitation temperatures. 



that from J = 2. This means that the ground and first ex- 
cited states of para-H2 can be thermahsed at lower densities 
compared to that of ortho-H2. In Fig. |1] we plot as a func- 
tion of temperature the critical hydrogen density for which 
the coUisional deexcitation rate becomes equal to the spon- 
taneous decay rate for the J=2— >J=0 transition. It is clear 
from this figure that the hydrogen density has to be high (in 
the range 60-175 cm"^) in order for the N{J=2)/N{3=0) 
ratio to be equal to the LTE value corresponding to typical 
kinetic temperatures inferred from the OPR (i.e 100 to 300 
K). 

In Fig.Elwe plot the ratio, N{J^2) / N{3^0), observed 
in DLAs, the Galaxy, LMC and SMC as a function of the 
total H2 column density. The vertical dotted lines in the fig- 
ure shows the expected values of the ratio for four different 
excitation temperatures assuming LTE. Values of the exci- 
tation temperature T02 for individual DLA H2 components 
obtained using Eq.|3are given in Table |21 The observed ex- 
citation temperatures are in the range 100 to 600 K with 
most of them at T02 — 150 — 300 K. If the level populations 
are in LTE then the required hydrogen density to maintain 
the equilibrium is 65 — 150 cm~'^ (see Fig. |1J. We can see 
from Fig. |S| that in DLAs where H2 is optically thick, the 
A'^(J=2)/A'^(J=0) ratio is larger than that seen in similar gas 
of the Galactic ISM, LMC and SMC. It can be seen from 
Table |5| that, in most of the DLAs, Toi is lower than or 



equal to T02 (see Table I^J. This is very much the case as 
well in most of the sightlines through the ISM and Magel- 
lanic clouds. It is well known that, due to a lower value of 
the Einstein coefficient of the J — 2 level compared to those 
of higher J levels, the UV and formation pumping processes 
can lead to enhancing the J = 2 level compared to the J = 
level. Thus the higher values of N{J=2) / N{J=0) seen in 
DLAs can be explained by higher pressure in the gas and/or 
higher radiation field. 

Fig. |n] gives the Ai'(J=3)/A''(J=l) ratio as a function 
of A''(H2). The vertical dotted lines in the figure shows the 
expected value of the ratio for four different excitation tem- 
peratures under the LTE assumption. The measurements in 
DLAs are consistent with local measurements and the ex- 
citation temperature T03 is in the range 100—680 K (see 
Table|21l. In Fig.|7|we plot the iV(J=3)/iV(J=l) ratio versus 
the A''(J=2)/A''(J=0) ratio. If formation and UV pumping 
contribute appreciably to populate the J = 2 and J = 3 
levels then we expect a tight relationship between the two 
quantities. The dotted line in the figure gives the expected 
relationship between the ratios under LTE. In the case of 
sightlines through the Galactic ISM, the LMC or SMC, the 
A'^(J=3)/A'^(J=1) ratio is higher than what is expected from 
the N{J=2)/N{J=0) value under LTE (or, T13 is higher than 
To2)- In the case of DLAs, most of the components have Tis 
close to T02 (points are on top of the dotted line). Note 
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Figure 6. The Af(J=3)/A'^(J=l) ratio observed in individual 
DLA components (black dots) is plotted against the total molec- 
ular hydrogen column density. The box drawn with dotted lines 
gives the allowed range for the 2abs= 1-973 components toward 
Q 0013— 004. Other data points are from Savage et al. (1977), 
Spitzer, Cochran & Hirshfelf (1974) for the Galactic ISM (tri- 
angles), and Tumlinson et al. (2002) for the LMC (squares) and 
SMC (asterisks) . The vertical dotted lines are the expected values 
of the column density ratio assuming the LTE for four different 
excitation temperatures. 

that these excitation temperatures are different from Tqi. 
This clearly means that UV pumping and formation pump- 
ing are not negligible even for the excitation of the low J 
levels. The nature of the local radiation field can be probed 
using excitations of Jjs 3 levels. This is what we do in the 
following Section. 

3.2.2 UV radiation field: High- J excitation 

It is known that in the photodissociation regions (PDRs) 
the J = 4 and J = 5 rotational levels are populated pre- 
dominantly by cascades following the formation of excited 
molecules and UV pumping from the low-J states. As ra- 
diative decay time-scales for these levels are very short com- 
pared to the coUisional time-scales, spontaneous decay is the 
main deexitation process. Among the two populating pro- 
cesses the UV pumping is an optical depth dependent pro- 
cess while the formation pumping is independent of optical 
depth. In an optically thick cloud, UV pumping is efficient 
in a thin shell surrounding the cloud. In the interior of the 
cloud, UV pumping becomes important only when the col- 
umn density becomes very large (i.e. through absorption in 
the damping wings). 

In Fig. |H| we plot log A'^(H2) as a function of 



Figure 7. The Ar(J=3)/Af (J=l) vs Af(J=2)/Af(J=0). DLA mea- 
surements are indicated by black dots. Other data points are from 
Savage et al. (1977), Spitzer, Cochran & Hirshfelf (1974) for the 
Galactic ISM (triangles), and Tumlinson et al. (2002) for the LMC 
(squares) and SMC (asterisks). The dotted line gives the expected 
relation under LTE with temperatures ranging from 100 to 600 K 
(horizontal tick-marks show the values for different temperatures 
with 50 K steps). 

log 7V(J=4)/7V(J=0) as measured in DLAs and along the 
ISM, LMC and SMC sightlines. As expected, a strong anti- 
correlation is present in the data, including DLAs. H2 ab- 
sorption lines in DLA systems have no strong overlapping 
wings. Therefore the high-J excitation is mostly due to 
photo- absorption in the systems with log A'^(H2) ^16.5 and 
to H2 formation in the systems with higher column densi- 
ties. Following analytic prescription by Jura (1975) we can 
write, 

P4,o/3(0)n(H2, J = 0)+ 0.2ARn{H)n = A(4 -> 2)n(H2, J = 4)(4) 

Here, /3(0), p4,o are, respectively, the photo-absorption rate 
in the Lyman and Werner bands and the pumping efficiency 
from J = to J = 4; A(4 2) is the spontaneous transition 
probability between J = 4 and J — 2 and R is the forma- 
tion rate of H2. Neglecting the second term in the left hand 
side of Eq. 2] leads to a conservative upper limit on the U V 
radiation field. The vertical dashed lines in Fig. |S] represent 
the corresponding predicted values of the A'^(J=4)/A'^(J=1) 
ratio for /3(0) = 2 x 10~'^°s~'^ (that is, approximatively the 
mean radiation field in the ISM) and f3{0) = 2 x lO'^'s-^ It 
can be seen from Fig.|H|that for log N{}i2) less than 16.5 the 
Ar(J=4)/A''(J=0) ratio in DLAs is of the order of or slightly 
higher than that seen in the ISM of our Galaxy. Quantita- 
tively the upper limits in most of the systems are consistent 
with 2 X 10"^° ^ P{0) < 2 X 10"^s"^ This probably means 
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Figure 8. The N{J=4:)/ N{J=0) ratio is plotted against the to- 
tal molecular hydrogen column density. DLA measurements arc 
indicated by black dots. Other data points are from Savage et al. 
(1977), Spitzer, Cochran & Hirshfelf (1974) for the Galactic ISM 
(triangles), and Tumlinson et al. (2002) for the LMC (squares) 
and SMC (asterisks). The box with dotted lines gives the allowed 
range for the ^abs= 1.973 component toward Q 0013— 004. The 
vertical dashed lines are the expected values of the ratio for two 
different values of the photo-absorption rate (/3o)- 



the optically thin H2 components without detectable H2 ab- 
sorption lines from the J = 4 state arise in gas embedded in 
a UV field with intensity similar to (or slightly higher than) 
that of the mean ISM field. 

There are two optically thin components in our sam- 
ple (2abs= 1.96822 toward Q 0013-004 and 3.02489 toward 
Q 0347—383) that show detectable J = 4 H2 absorption lines. 
Detailed analysis of these component suggests an ambient 
field intensity consistent with few times the mean ISM field 
intensity (Petitjean et al. 2002; Levshakov et al. 2002). The 
same conclusion was derived by Reimers et al. (2003) for 
the optically thin H2 component at z^^s— 1-15 system to- 
ward HE 0515 - 4414. 

The above ratio has similar values at high log A'^(H2) 
in DLAs and in our Galaxy. This is a hint for the forma- 
tion pumping in DLAs with high A''(H2) being similar to 
the local one. There are two optically thick components (at 
«abs= 2.59471 toward Q 0405-443 and Zabs= 2.08696 to- 
ward Q 1444-1-014) that do not show detectable absorption 
lines from the J = 4 state. In the latter system the ratio 
Af(J=4)/iV(J=0)^ 10"". This is much lower than the val- 
ues seen in the ISM at similar total A'^(H2) and could be a 
consequence of lower H2 formation rate in this system. High 
values of the radiation field intensity were inferred for some 
of the optically thick components when the contribution of 
the second term in Eq. is estimated using the average 
metallicity and dust depletion (Ge & Bechtold 1997; Petit- 
jean et al. 2000; Ge, Bechtold & Kulkarni 2001). 



4 ANALYSIS OF CARBON ABSORPTION 
LINES 

As the ionization potential of C I is 11.2 eV, the ionization 
state of Carbon is sensitive to the same photons that destroy 
H2. Therefore, C I is usually a good tracer of the physical 
conditions in the molecular gas (see however Srianand & Pe- 
titjean 1998). In what follows we investigate the relationship 
between the detectability of C I absorption line and other 
measurable quantities in our spectra. We derive additional 
constraints on the physical conditions in DLAs using C I 
fine-structure absorption lines. 



4.1 Detectability of C I absorption lines 

The results of simultaneous Voigt profile fitting to C I, C I* 
and C I** absorption lines in our sample are summarised in 
Table 01 

In the interstellar medium of our Galaxy, all clouds with 
log N{}i l) ^ 21 have log iV(H2) > 19 and log iV(C l) > 14 
(Jenkins & Shaya 1979; Jenkins, Jura & Loewenstein 1983). 
In our sample, C I absorption lines are detected in most of 
the DLAs that show H2 absorption lines(see also Ge & Bech- 
told 1999). There are three exceptions: the components at 
Zabs= 2.59471 and 2.59486 toward Q 0405-443 and at Zabs = 
2.81100 toward PKS 0528-250. These C l non-detections are 
surprising as the H2 absorption lines from these components 
are optically thick so that C I is expected to be conspicuous. 

Usually, DLAs in which no H2 is detected through the 
whole profile do not show any detectable C I absorption line 
(with a typical upper limit of 10^^ cm~'^). The only excep- 
tion is the high-metallicity sub-DLA at 2:abs= 2.139 toward 
Tol 1037-270 (see Srianand & Petitjean 2001). On the con- 
trary, in DLAs where H2 is detected, some components show 
detectable C l absorption line without detectable molec- 
ular absorption (A'^(H2)^ 10^*cm~^). This is the case in 
Q 0013-004 (Petitjean et al. 2002) and Q 0551-366 (Ledoux 
et al. 2002). 

Note that C i is also detected at 2abs= 2.28749 toward 
Q 2332—094 but the presence of H2 molecules can not be 
probed in this system due to the presence of an interven- 
ing Lyman limit system. The sub-DLA at 2; = 1.15 toward 
HE 0515—4414 shows C I and H2 absorption lines (Quast et 
al. 2002; Reimers et al. 2003). C I absorption lines have also 
been detected at 2abs= 1.776 system toward Q 1331-1-170 
(Chaffee et al. 1988). Presence of H2 in this system is re- 
cently reported (Cui et al., 2004). 



4-1.1 Dependence on H2 column density 

Jenkins & Shaya (1979) found 7V(C l) does not scale linearly 
with either of Af(H l), iV(H2) or iV(Htotai) in the Galac- 
tic ISM. They explained this behavior as a result of strong 
differences in the response of C I, H I and H2 to physi- 
cal conditions (electron density, temperature etc.), coupled 
with marked variations of these conditions from one cloud 
to the other. In Fig. |U] we plot the C I column density as 
a function of H2 column density in individual components. 
Among the systems that show H2 absorption lines(filled cir- 
cles with error-bars) there is no clear trend between A''(H2) 
and A'^(C l) even though the presence of C l absorption lines 
usually indicate the presence of H2 (see discussion above) . 
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Table 3. DLA systems with detected Cl absorption lines in the sample of Ledoux et al. (2003) 



Quasar/ 
log N(H I) 






logiV(/f2) 


log iV(Cl) 


logAf(Ci*) 


logiV(Cl**) 


fa(km s~l) 


Tcmb(i^) 


nH(c 


n,-3) 


P/k (cm~3 K) 




Q 0013-004 
20.83(0.05) 


1 
1 


96679 
96691 




13.06 ± 0.04 
12.86 ± 0.03 


12.72 ± 0.05 
12.55 ± 0.07 




3.1±0.6 
3.0±0.2 


< 
< 


13 
14 


8 


20 
20 


- 70 

- SO 


2000-7000 
2000-8000 




1 


96706 




12.31 ± 0.07 


12.26 ± 0.11 




4.93±1.44 
















1 


96822 


13.04 ± 0.01 


13.10 ± 0.01 


12.77 ± 0.02 


4.3±0.2 


< 


20 





170 


- 200 


11730 - 12400 




1 


97280 


''17.48-19.78 


13.08 ± 0.01 


12.77 ± 0.03 




6.0±0.3 


< 


13 





40 


- 60 


4000-6000 




1 


97296 




13.44 ± 0.01 


13.20 ± 0.01 


12.53 ± 0.04 


6.2±0.2 


< 


13 




50 


- 65 


5000 — 6500 










-Li.OD It U.llO 


i-A.OO It U.±U 
















±UUU — OOl7l7 




1 


97365 


'^17.30—19.60 


12.32 ih 0.05 


< 12.11 




5.7+0.4 




15 







85 


< 8500 




1 


97382 




12.40 ± 0.04 


12 .04 lb 0.13 




4.5 + 0.4 




15 






95 


< 9500 




1 


97399 




11.91 ± 0.13 






7.3+1.0 
















1 


97417 




11.49 ± 0.25 






0.4 + 0.3 
















1 


96737 


< 14.0 


12.91 ± 0.03 


12.77 ± 0.06 




14.4+1.1 


< 


17 




40 - 


- 135 


4000-13500 










±^.oo in u.uo 


12 68 i 09 








19 


3 










1 


97109 


< 14.0 


12.28 ± 0.03 


< 12.11 




10.0 + 0.9 


< 


21 





< 


200 


< 20000 




1 


97144 


< 14.0 


12.96 ± 0.03 


12.91 ± 0.03 




13.4 + 0.9 


< 


18 




65 - 


- 160 


6500-16000 


Q 0347-383 
20.73(0.05) 


3 


02485 




11.73 ± 0.26 


< 11.50 


< 11.75 


4.9 + 0.3 


< 


17 





< 


48 


< 4800 


Q 0405-443 


2 


59474 


,,,,+0.07 


< 12.23 




















21.05(0.10) 


2 


59485 


15 Sl + O-lS 


< 11.90 




















Q 0528-250 


2 


81100 


17.93 + °-^„ 


< 12.00 




















21.35(0.10) 


2 


81112 


17.90+0-;^ 
< 14.0 


12.36 ± 0.10 


12.30 ± 0.10 




0.6 + 0.1 


< 


17 





25 - 


- 270 


3250 - 17000 


Q 0551-366 


1 


96152 


12.69 ± 0.07 


< 12.18 


< 11.94 


4.3 + 1.4 


< 


10 


5 


< 


17 


< 1700 


20.70(0.10) 


1 


96168 


15.80 + "-- 
< 14.0 


12.64 ± 0.07 


12.84 ± 0.07 


< 12.16 


2.1 + 0.7 


< 


22 


3 


170 


- 185 


12950 - 14280 




1 


96180 


12.42 ± 0.13 


< 12.18 


< 11.94 


3.9 + 2.3 


< 


12 


9 


< 


56 


< 5600 




1 


96214 


"■«tS:r3 

i5.asl°- 

< 14.0 


12.66 ± 0.12 


12.69 ± 0.11 


12.11 ± 0.34 


2.1 + 0.8 


< 


18 


7 


55 - 


- 390 


8250 - 30400 




1 


96221 


13.16 ± 0.06 


12.98 ± 0.09 


12.26 ± 0.36 


12.8 +1.7 


< 


14 




30 - 


- 150 


3840 - 15150 




1 


96268 


12.63 ± O.OS 


< 12.18 


< 11.94 


4.0 + 1.8 


< 


10 





< 


25 


< 2500 


Q 1037-270 


2 


13900 


< 14.0 


12.51 ± 0.02 


< 12.48 




15.5 + 0.6 


< 


16 


4 


< 


105 


< 10500 


19.70(0.05) 


2 


13940 


< 14.0 


12.48 ± 0.02 


< 12.30 




4.5+1.2 


< 


13 


2 


< 


65 


< 6500 


Q 1232+082 


2 


33771 


19.57+0- 


13.86 ± 0.22 


13.43 ± 0.07 


12.63 ± 0.22 


1.7 + 0.1 


< 


14 


40 


- 60 


2320 - 5280 


20.90(0.08) 


























Q 1444+014 


2 


08679 


16.60+°'^'* 


12.67 ± 0.09 


12.52 ± 0.14 




11.6 + 2.8 


< 


18 




20 - 


- 110 


6560 - 27280 


20.25(0.07) 


2 


08692 




12.82 ± 0.11 


12.42 ± 0.12 




1.1 + 0.3 


< 


14 





4 - 


- 54 


824 - 10260 


Q 2332-094 


2 


28749 


13.33 ± 0.03 


13.14 ± 0.03 


12.34 ± 0.11 


3.8 + 0.3 


< 


14 


6 


45 


- 75 


4500 - 7500 



20.25(0.07) 

" total H2 column density in tliree components at 2j^bs= 1-96679, 1.96691 and 1.96706 is quoted, 
total H2 column density in three components at 2abs= 1 97280, 1.97296 and 1.97316. 

total H2 column density in four components at aabs= 1 97365, 1.97382, 1.97399 and 1.97417 is quoted. 



4-1.2 The Carbon ionization state 

The probability of detecting C I is expected to be higher 
in systems with higher A''(H l) and/or metallicity. Ideally, 
we would like therefore to know A''(H l) for each individual 
C I components. This is not possible as all components are 
definitely blended in one strong H I DLA absorption line. 
Estimation of N(H l) is possible when the H2 component 
is well separated from the rest of the components (as in 
2abs= 1.96822 toward Q 0013-004) or when 21 cm obser- 
vations are available (as in the case of H2 components to- 
ward Q 0528—250). Note that the presence of very strong C I 
absorption line in the component at Zabs^ 1.96822 toward 
Q 0013-004 (that has log N(H l)s;i9.4) is mainly due to high 
metallicity (Petitjean et al. 2002). Whereas the absence of 
C I in the component at 2abs= 2.81100 and the weakness of 
C I line of the component at 2abs= 2.81112 component to- 
ward Q 0528—250 are probably due to excess radiation field 
from the QSO. 

The ionization state of Carbon is difficult to determine 
as the C nA1334 absorption line is usually highly saturated. 
We can however partly overcome this difficulty assuming 
that conditions are fairly homogeneous in the DLA system. 
Under the assumptions that the enrichment of Carbon fol- 
lows that of Of- elements and that the relative depletion be- 
tween Sulfur and Carbon is negligible, we can use the well 
determined A''(S 11) as an indicator of A'" (C 11) (see Fig. [101. 
In the sun, the Carbon abundance is 1.28 dex higher than 
that of Sulfur and typical depletion of Carbon relative to 
Sulfur in the Cold ISM is 0.4 dex (see Table 5 of Welty et 



al. 1999). The two dotted lines in Fig. 1101 give the expected 
correlation for log N{C l)/N{C 11) = —3 (lower line) and —2 
(upper line) respectively, when relative solar abundances are 
used and it is assumed that there is no depletion of Carbon 
relative to Sulfur. The short-dashed lines give the same cor- 
relations when a depletion of Carbon relative to Sulfur of 0.4 
dex is further assumed. If the absorbing gas originates from 
the CNM then log iV(C i)/A^(C 11) is expected to be more 
than —3 (see Fig. 3 of Liszt 2002). Therefore, within uncer- 
tainties due to depletion, it is apparent from Fig. 1101 that 
the DLA components with C I detections have a ionization 
state consistent with them originating from the CNM. 

It is interesting to note that the distribution of A'^(S 11) 
is somewhat similar for components with both H2 and C I 
absorption lines (filled circles), and for components with C I 
but no II2 absorption lines (open circles). However C I col- 
umn densities are typically lower in the components with- 
out H2 suggesting, as expected, that when H2 is seen, the 
C l/C II ratio is larger. 

Most of the upper limits on C I are consistent with 
iV(C i)/iV(C 11) ^ -3 (see Fig. Unj. This can mean either 
that the relative depletion of Carbon compared to Sulfur is 
larger than 0.4 dex in the CNM, which is unlikely, or that 
most of the DLA systems originate from the warm neutral 
medium (WNM) or warm ionized medium (WIM) where the 
above ratio can be as low as 10"*. 
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Figure 9. The column density of C I measured in individual 
components of DLAs is plotted as a function of N(il2). The filled 
circles with error-bars are points from the components that show 
detectable H2 absorption lines. The open circles are the measure- 
ments from components that show C I without detectable H2. 



4.1.3 The effect of dust 

In Fig. 111! we plot log A'^(C l) against the depletion fac- 
tor defined as log (iV(Fe ll)/iV(X n))-[Fe/X]0 with either 
X = Zn, S or Si. C l absorption line is not detected in sys- 
tems with low depletion factors (i.e., [Fe/Zn] less than —0.5 
dex) whereas components with higher depletion factors read- 
ily show detectable C l absorption lines. This trend is not 
surprising as there is a 4a correlation between the deple- 
tion factor and the metallicity of the gas in our sample (see 
Figure 12 of Ledoux et al. 2003). Depletion factors lower 
than 0.5 dex are usually seen in systems with [Zn/H] ^ —1. 
High depletion factor in high metallicity gas implies high 
dust content and hence high dust optical depth to the UV 
radiation. The absence of C I in components with low dust 
depletion is a combination of low metallicity and low dust 
optical depth to the UV radiation. It is worth remembering 
that similar relation exists between the detectability of H2 
and depletion (Fig. 14 of Ledoux et al. 2003). 



4.2 C I fine-structure excitation 

In most of the DLAs with C I detections we also detect 
absorption lines from the excited fine-structure levels. It is 
therefore possible to use the relative populations of the C I 
ground state levels to discuss the particle density, the am- 
bient UV radiation field and the temperature of the cosmic- 
microwave background radiation (see Bahcall et al., 1973; 
Meyer et al., 1986; Songaila et al. 1994; Ge, Bechtold & 
Black, 1997; Roth & Bauer, 1999; Srianand et al. 2000; Silva 
& Viegas, 2002; Quast et al. 2002). In the Galactic ISM, 
fine-structure excitation of C I has been used to study the 
distribution of thermal pressure (see Jenkins & Tripp 2001). 
In Fig.IUlwe plot the ratio N{C l*)/Af(C l) as a func- 



14 



" 13 



12 



1 1 



' 1 ' 




1 








_ __.'r^ 










- 









13 



13.5 



14 14.5 
log N(S 11) (cm-2) 



15 



Figure 10. The C I column density measured in individual com- 
ponents of DLAs is plotted as a function of the S II column den- 
sity. The filled circles with error-bars are points from components 
that show detectable H2 absorption lines. The open circles are the 
measurements from components that show C I but no detectable 
H2. The arrows are the upper limits from systems where there 
is no H2 and/or C I absorption lines. The two dotted lines give 
the expected correlation for log N(C l)/N(C 11) = —3 (lower line) 
and —2 (upper line) respectively, when solar relative abundances 
are used and it is assumed that there is no depletion of Carbon 
relative to Sulfur. The short-dashed lines give the same correla- 
tions when a depletion of Carbon relative to Sulfur of 0.4 dex is 
further assumed. The C I detections are consistent with what one 
expects in the case of the CNM. 



tion of A'^(C l). It is clear that the C I column densities in 
DLAs (filled circles with error-bars) are at least an order 
of magnitude less than that measured in the ISM (stars). 
This is probably a consequence of lower metallicities and/or 
low H I content in DLA components. The important point is 
that the N{C l*)/N{C l) ratio measured in DLAs is remark- 
ably larger than in the Galaxy. However, while comparing 
the ISM and DLAs, it is important to remember that most 
of the sightlines used by Jenkins & Tripp (2001) have H2 
fraction orders of magnitude higher than what we measure 
in DLA components. As H2 collisions are less efficient in 
populating the excited fine-structure state of C I we expect 
that for a given total hydrogen density (and a given kinetic 
temperature) N{C l*)/N{C l) be higher in DLAs. 

The horizontal dotted line in Fig. 1121 indicates the ex- 
pected value of the ratio if it is assumed that the excitation 
is due to the CMBR only with a temperature Tcmbr = 8.1 
K as expected at 2 = 2, the typical redshift of our sample. 
It is clear that the CMBR field expected from the Big-Bang 
is not sufficient to explain the observed ratios and an extra 
contribution is required from coUisional processes and/or the 
UV flux. 

Collisions with H, He, e, p and H2 can populate the 
excited state of C I. As H I is the dominant form of hy- 
drogen in the gas, the contribution to the fine-structure 
excitation by H2 collisions can be neglected. The electron 
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Table 4. Physical conditions in systems in which C I is detected 



QSO 






log Af(S II) 


log N{C II*) 


ric (cm •^)" 


njj (cm 3) 


Q 0013 — 


004 


1.96679+ 


14.27 ± 0.05 




2.7 10~ 








1.96822^^ 


14.16 ± 0.02 




4.9 10"'^ 








1.97280+ 


14.15 ± 0.06 




3.8 10~^ 








1.97296^ 


14.75 ± 0.03 




2.3 10~^ 








1.97316+ 


13.91 ± 0.05 




2.4 10~^ 








1.97365+ 


13.86 ± 0.08 




1.5 10~^ 








1.97382^ 


14.12 ± 0.05 




8.0 10~^ 








1.97399+ 


13.55 ± 0.16 




6.9 10~^ 








1.96737 


14.24 ± 0.05 




2.4 10~^ 








1.96763 


14.12 ± 0.07 




2.7 10"'^ 








1.97109 


14.77 ± 0.01 




1.0 10~^ 








1.97144 


14.08 ± 0.03 




4.3 10~^ 




Q 0347 — 


383 


3.02463 


14.50 ± 0.03 


< 13.00 


< 5.0 10"* 


< 6.3 






3.02485+ 


14.36 ± 0.04 


13.55 ± 0.23 


2.3 10"'' 


4.4 — 41.9 






3.02501 


13.51 ± 0.12 


< 12.72 


< 4.9 10"'^ 


< 32 


Q 0405 — 


443 


2.59440 


13.91 ± 0.04 


< 12.41 


< 3.7 10~^ 


< 6.3 






2.59464 


14.72 ± 0.02 


13.09 ± 0.10 


< 5.7 10"* 


3.0 — 7.5 






2.59474+ 


14.07 ±0.09 


13.20 ±0.20 


< 2.5 10-3 


7.7 - 57.3 






2.59485+ 


14.73 ±0.02 


13.25 ±0.08 


< 5.6 10"* 


4.5 - 10.3 


Q 0551 - 


366 


1.96152 


14.60 ±0.10 


< 12.00 


3.7 10-3 


< 0.5 






1.96168+ 


14.01 ±0.27 


13.93 ±0.09 


3.2 10-2 


123.0 - 331.8 






1.96180 


14.23 ±0.19 


< 12.00 


4.6 10-3 


< 1.2 






1.96214+ 


Blended 


12.58 ±0.09 










1.96221+ 


Blended 


13.36 ±0.07 










1.96268 


14.63 ±0.11 


12.98 ±0.08 


3.0 10-3 


2.4-9.6 


Q 1037 - 


270 


2.13900 


14.68 ±0.02 


< 13.50 


2.0 10-3 


< 13.2 






2.13940 


14.23 ±0.11 


13.31 ±0.01 


5.3 10-3 


< 14.2 


Q 1232 + 


082 


2.33771+ 


15.24 ±0.11* 


< 14.00 


< 7.0 10-2 


< 48 


Q 1444 + 


014 


2.08679+ 


14.32 ±0.05 


13.12 ±0.08 


1.1 10-2 


0.9 - 9.9 






2.08692+ 


14.17 ±0.06 


12.78 ±0.20 


1.9 10-2 


< 13.0 



A label "+" indicates a component in which H2 is detected 
° Derived from the Carbon ionization equilibrium 
column density of Si II is given. 



and proton densities are expected to be very small, at least 
smaller by two orders of magnitude than the hydrogen den- 
sity, and their contribution is also negligible (Keenan et al. 
1986). The He I coUisional rates are much less than that of 
H I and the He l/H l ratio is small which makes collisions 
with He I unimportant (see Fig. 1 of Silva & Viegas 2002). 
Thus in our analysis of the C I excitation we consider only 
collisions by neutral hydrogen. The rates are taken from 
Launay & Roueff (1977). The spontaneous decay rates are 
Aio = 7.93 X 10"** s"^ and A21 = 2.68 x 10"'^ s"^ (Bahcall 
& Wolf 1968). The corresponding CMBR excitation rate is 
derived from these values. The UV pumping rate in the cloud 
depends on the nature and strength of the UV radiation 
field. We assume that the UV intensity is the same as in the 
ISM of our Galaxy as suggested by the high— J excitation 
of H2. 

The dashed lines in Fig. [H give the expected ratio for, 
TcMBR = 8 K, a UV radiation field like in the Galaxy (with 
an excitation rate of Foi = 7.55 x IQ-^" s-^), hydrogen 
density in the range nu ~ 20—250 cm-3, versus the kinetic 
temperature, Tkin, in the range 80 — 200 K (consistent with 
Toi measured in DLAs). It can be seen that a typical den- 
sity range consistent with most of the observed point is, 
20< nH(cm-3) ^150. The density (nn) and pressure {p/K) 
derived for individual components are summarized in col- 
umn 9 and 10 respectively in Table. |3 Here, we assume 



TcMBR = 2.7x(l+Zabs) and T — Toi in the case of H2 de- 
tection and T — 100 K otherwise. The derived pressure range 
in DLA components are higher than that typically measured 
in the galactic ISM (see Jenkins & Tripp 2001) and consis- 
tent with what is expected in a the cold neutral medium 
(CNM) with lower metallicity (Z~0.1 Z©) and dust deple- 
tion (see Wolfire et al. (1995, 2001), Wolfe et al., (2003a, b) 
Srianand et al. 2005) . Interestingly the derived density range 
in most of the components is close to the critical density for 
thermalising the H2 iV(J=2)/iV(J=0) ratio (see Fig.gJ. 



5 CARBON IONIZATION STATE 

In the ISM, when hydrogen is neutral, most of the carbon 
is in the form of C 11 so that we can use the C l/C 11 ra- 
tio to derive the physical conditions in the gas. As already 
mentioned, the C 11 absorption features are always badly 
saturated and we rely on the realistic assumption that C 11 
in the neutral phase can be traced by S 11 and/or Si 11 so 
that we can use the weak lines of these species to derive the 
C II column densities in the components of interest. 

Assuming photoionization equilibrium between C I and 
C II and using the atomic data from ShuU & Van Steenberg 
(1982) and Pequignot et al. (1991), we can write. 
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-0.5 -1 -1.5 -a 

[Fe/7n] or [Fe/Si] or [Fe/S] 

Figure 11. The column density of C I measured in individual 
components is plotted as a function of the metal depletion fac- 
tor computed using [Fe/Zn] or [Fe/S] or [Fe/Si]. It is apparent 
from this figure that C I absorption line is usually detected in 
components with high depletion factors. 
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Figure 12. The ratio A'"{C l*)/N{C l) is plotted on a logarithmic 
scale as a function of A'^(C l). The points with error bars are 
our measurements in individual DLA components. The stars are 
data measurements in the ISM of the Galaxy drawn from Jenkins 
& Tripp (2001) and Jenkins, Jura & Loewenstein (1983). The 
horizontal dotted line gives the expected value of the ratio if it 
is assumed that C I is excited by the CMBR only with Tcmbr 
= 8.1 K as expected at z = 2. The short dashed lines give the 
expected ratio for different nj] as a function of temperature (top 
portion of the x-axis). 




log (Electron density (cm"'*)) 

Figure 13. Histogram of the electron density derived using Eg. 151 
for T = 100 K. The shaded histogram is based on systems that 
show detectable C I absorption lines. The histograms with con- 
tinuous and dotted lines are upper limits on the electron density 
for, respectively, systems with log 7V(C ll) greater and less than 
15. We use the C I photoionization rate given for the Galactic UV 
field by Pequignot & Aldrovandi (1986). 



(5) 



r iV(CII)l,104 
Here, F is the photoionization rate for C I. In the local ISM, 
Tgai ~ 2 - 3.3 X 10~^° s'^ (Pequignot & Aldrovandi 1986). 
Here, we neglect the ion-molecular interaction and charge 
exchange reactions that may produce C I. Thus rie/r can 
be constrained once the temperature of the gas is known. 
Application of Eq. to the cold Galactic ISM (with T = 
100 K) has resulted in rio ~ 0.14 ± 0.07 cm"^ (Welty et al. 
2002 and references there in). In the stable CNM considered 
by Wolfire et al. (1995) the electron density is expected to 
be in the range 0.01^ nc(cm"^) <0.02 for Z = 0.1 Zq and 
dust abundance one tenth of the Galactic ISM. 



5.1 Systems with C I detections 

First we concentrate on the systems with C I detections. We 
estimate the electron density assuming a UV field similar 
to the Galactic mean field (i.e., Fgai = 2.5 x 10"^^° 
T = Toi for the H2 components and T =100 K otherwise. 
Individual values of rio derived for these systems are given 
in Table 2] The electron density is in the range 0.7x10-2 ^ 
rio (cm"^) s; 4.9 x lO^^. Together with uh given in TableOl 
this suggests that ric/nn ^ 10"^ for most of the systems. 
Therefore, the ionization state of the gas with H2 and C I is 
similar to that in the CNM in a moderate radiation field. 



5.2 Systems without H2 and C I 

In the case of systems in which neither H2 nor C I are de- 
tected, we assume T = 100 K and Fg^i = 2.5 x 10"^° s"^ 
and obtain upper limits on the electron density. The results 
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are plotted in Fig. 1131 We notice that the inferred electron 
densities are much smaller than in systems in which C I 
(and H2) are detected. The difference is even larger if we 
use only the systems with log A'^(C 11) ^15. For this gas the 
inferred electron densities are less than 10~^ cm^'' (with a 
median of 10"'' cm"''). This may indicate that the absorp- 
tion originates from warm neutral medium (say T = 8000 
K). If the average radiation field in DLAs is similar to that 
in the Galactic ISM then the absence of C I in most of the 
DLAs could just be a consequence of lower densities (and/or 
higher T) in these systems. One can derive an independent 
estimate of the particle density using the excitation of C 11 
fine structure levels. This is what we do in the following 
Section. 



6 EXCITED FINE-STRUCTURE LINE OF C II 

6.1 Method to derive physical parameters 

Under LTE, the column density ratio 7V(C ll*)/N{C 11), 
can be written as, 

iV(C ir) ^ Qi2(e) ue + Qi2(H) uh + ri2(CMB) 
N{C II) A21 ^ ' 

where, Qi2(e) = 7.8 x 10-^exp[-91.27/r] T-"-^ cm"^ s'^ 
and Qi2(H) = 1.3 x 10~^ea;p[-91.27/r] cm'^* s"^ are the 
collisional excitation rates per unit volume for electrons and 
hydrogen atoms (Bahcall & Wolf 1968) with T being the 
kinetic temperature of the gas. The Einstein's coefficient is 
.421 = 2.291x10"'' s-\ The CMB pumping rate, ri2(CMB), 
equals 6.6x10"^^ s~^ and 1.1x10"^ s~^ for, respectively, 
redshifts 2 and 3 (Silva & Viegas 2002). For a given temper- 
ature, the collisional excitation rate for electrons is orders 
of magnitude larger than that for hydrogen atoms. For ex- 
ample, when T — 1000 K, whenever n^/nu is larger than 
5xl0~^, collisions with electrons is the dominant process. 
UV pumping is an additional possible excitation mechanism. 
For the mean radiation field in our Galaxy, the UV pumping 
rate is 9.3x10"" s"^ (Silva & Viegas 2002). This is simi- 
lar to or slightly lower than the CMB pumping rate for the 
range of redshift we consider in this study. 

We compute the expected value of log N{C ii*)/N(C 11) 
ratio under different situations and in particular the WNM 
and CNM solutions given in Table 3 of Wolfire et al. (1995). 
Note that if CMBR pumping alone is responsible for the 
excitation, the expected ratios are —4.54 and —3.31 for, 
respectively, Zabs= 2 and 3. For the standard ISM (with 
stable pressure in the range 990—3600 cm~^ K), we de- 
rive -2.62 ^log N{C n*)/iV(C 11)^ -2.20 for the CNM 
and -3.65 ^log iV(C ll*)/N{C 11)^ -3.17 for the WNM. 
As the DLA gas has low metallicity and low dust con- 
tent, the expected pressure should be higher in the two 
DLA phases (see Liszt 2002 and Wolfe et al. 2003). If 
we assume Z — O.IZq and dust to gas ratio one tenth 
of the ISM value (which is typical of DLAs) then we ex- 
pect -2.26 sClog iV(C n*)/7V(C ii)s; -2.08 for the CNM 
and -3.39 <log N{C n*)/iV(C 11)^ -2.70 for the WNM 
in DLAs. Thus if DLAs originate from H I gas in a two- 
phase equilibrium, we expect the CMBR pumping to be 
sub-dominant compared to collisional excitation. In fact if 
the gas is completely neutral then from Eqs. (jSJ and JSJ we 



derive, 

A2lN(C IF) aeUc 

(7hN(C II) an 
A2iN(C IF) _ acN(C I)r 
aHN(C II) crHarN(C II) ' ^ ' 

This gives an independent estimate of uh which, by com- 
parison with the estimate derived from the C I excitation, 
can lead to constraints on the radiation field. However, the 
ionization fraction (i.e ria/nu) and the temperature of the 
neutral absorbing gas must be accurately determined be- 
fore densities can be derived using the A'^(C ll*)/A'^(C 11) 
ratio. Indeed, the fact that, in DLAs, the Al ill absorption 
profile is very similar to that of neutral or singly ionized 
species has been used as evidence for the presence of ion- 
ized gas being mixed with the neutral gas in DLAs (Lu et 
al. 1996; Prochaska & Wolfe 1999; Howk & Sembach 1999; 
Wolfe & Prochaska 2000; Vladilo et al. 2001; Izotov et al. 
2001). In the ionized gas (i) well as T will be higher 

than what is expected in the warm or cold neutral gas and 
(ii) A'^(Si 11) will under-predict A'^(C 11) as the ionization cor- 
rections are different for the two species (see Fig. 1 of Izotov 
et al. 2001). Neglecting the presence of ionized gas can arti- 
ficially enhance the derived uh values. 

6.2 Frequency of C 11* detection 

In our sample, all the systems that show H2 also show de- 
tectable C II* absorption line. The details of the fits to 
the G II* absorption line for these systems are summarized 
in Table |S] and shown in Fig. 1141 The components toward 
Q 0013 - 004 and Q 0528 - 250 are badly blended and it is 
therefore not possible to fit A''(G 11*). We detect G 11* in 
seven out of the 21 DLAs that do not show II2 absorption 
lines. If we also include the 8 DLAs that show H2, about 50% 
(15 out of 29) of the DLAs in our sample show detectable 
G II* absorption line which is consistent with the finding by 
Wolfe et al. (2003). 

In Fig. 1151 we plot the average (over the whole profile) 
N{C ll*)/N{C 11) ratio measured in DLAs in our sample to- 
gether with the measurements by Wolfe et al. (2003a, 2003b) 
against the total H I column density and the silicon metal- 
licity. Here we use the total column density summed over all 
the components. A'^(C 11) is computed from A'^(Si 11) assum- 
ing solar abundance ratio without any ionization correction. 
From the upper panel in Fig. 1151 it can be seen that G 11* 
is detected in all the systems with log A'^(H l) ^21.0. In- 
terestingly no such relationship exists between C I (or H2) 
and H I. Most the systems with log N{R l) ^21.0 have 
N{C ll*)/N{C 11) consistent with what is expected in CNM. 
On the contrary, the measured values of N{C ll*)/N{C 11) in 
systems with lower A''(H l) spread over more than an order 
of magnitude covering the expected ranges for WNM and 
CNM. 

From the bottom panel, it can be seen that G 11* is 
frequently detected in gas with high metallicity as already 
noticed by Wolfe et al. (2003a). Most of the systems that 
show G II* absorption line with lower N{}i l) do have statis- 
tically higher metallicity. In the whole sample the number 
of systems with G 11* detections that are consistent with 
CNM and WNM are approximately equal. Most of the up- 
per limits on the ratio, measured in the metallicity range 
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Table 5. Results of Voigt profile fits to DLAs that do not show detectable H2 absorption lines 



Z/Zq 



logJV(Cl) 



[Fc/X] 



log iV(C I] *) 



Q 0000- 
Q 0010- 


-263 
-002 


21.40(0.08) 
20.95(0.10) 


-2.06(0.09) 
-1.43(0.11) 


Q0058- 


-292 


21.10(0.10) 


-1.53(0.10) 


Q0102- 


-190 


21.00(0.08) 


-1.90(0.09) 


Q0112- 


-306 


20.50(0.08) 


-2.43(0.09) 



20.30(0.10) 



Q 0112-1-029 20.90(0.10) 



Q0135-273 21.00(0.10) 



Q0405-443 21.15(0.15) 



20.45(0.10) 



Q0841-I-129 21.05(0.10) 
20.80(0.10) 

Q 1037-270 19.70(0.05) 
Q 1101-264 19.50(0.04) 



Q 1117-134 20.95(0.10) 
Q 11574-014 21.80(0.10) 



Q 12234-178 21.40(0.10) 

Q 13374-113 20.12(0.05) 
21.00(0.08) 

Q 14514-123 20.40(0.10) 
20.20(0.20) 

Q 2059-360 20.29(0.07) 
20.98(0.08) 



Q 2138-444 20.98(0.05) 
Q 2332-094 20.50(0.07) 



-0.50(0.15) 



-1.32(0.15) 



-1.40(0.10) 



-1.41(0.15) 



-2.04(0.11) 



-1.59(0.10) 
-1.60(0.10) 

-0.31(0.06) 
-1.07(0.06) 



-1.42(0.14) 
-1.40(0.10) 



-1.63(0.11) 

-1.81(0.07) 
-1.86(0.10) 
-2.27(0.14) 
-2.10(0.21) 

-1.94(0.08) 
-1.78(0.11) 



-1.74(0.05) 
-1.33(0.11) 



3.39013 
2.02471 
2.02484 
2.67123 
2.67142 
2.36958 
2.36966 
2.41844 
2.41861 
2.41869 
2.70111 
2.70163 
2.70217 
2.70233 
2.70257 
2.70276 
2.70316 
2.70332 
2.70334 
2.42234 
2.42248 
2.42266 
2.42277 
2.42298 
2.42323 
2.42338 
2.42353 
2.79948 
2.79960 
2.79972 
2.79982 
2.79998 
2.80010 
2.80025 
2.54967 
2.54973 
2.54992 
2.55011 
2.55049 
2.55060 
2.55078 
2.55087 
2.55097 
2.55105 
2.55124 
2.62101 
2.62140 
2.62158 
2.62178 
2.62219 
2.62235 
2.62287 
2.62311 
2.62320 
2.37439 
2.37455 
2.47604 
2.47621 
2.47642 
2.13910 
2.13951 
1.83817 
1.83831 
1.83854 
1.83871 
1.83890 
1.83915 
1.83925 
1.83933 
3.35027 
3.35046 
3.35067 
1.94311 
1.94318 
1.94354 
1.94346 
1.94361 
1.94375 
1.94403 
1.94385 
1.94426 
2.46530 
2.46559 
2.46607 
2.46628 
2.50766 
2.50792 
2.79567 
2.79584 
2.46897 
2.46921 
3.17081 
3.17112 
3.17142 
2.50734 
2.50753 
3.08261 
3.08291 
3.08314 
2.85234 
3.05632 
3.05657 
3.05676 
3.05690 
3.05722 
3.05738 



< 12.15 

< 12.15 

< 12.35 

< 12.35 

< 12.30 

< 12.30 

< 12.32 

< 12.32 

< 12.32 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 12.20 

< 11.78 

< 11.78 

< 11.78 

< 11.78 

< 11.78 

< 11.78 

< 11.78 

< 12.22 
<12.00 
<12.00 
<12.00 
<12.00 
<12.00 
<12.00 
<12.00 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.70 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 11.64 

< 12.21 

< 12.21 

< 12.21 

< 12.21 

< 12.21 
12.51(0.02) 
12.46(0.02) 

< 11.68 

< 11.68 

< 11.68 

< 11.68 

< 11.68 

< 11.68 

< 11.68 

< 11.68 

< 12.76 

< 12.45 

< 12.45 

< 12.15 

< 12.15 

< 12.15 

< 12.15 

< 12,15 

< 12.15 

< 12.15 

< 12.15 

< 12.15 

< 12.30 

< 12.30 

< 12.30 

< 12.30 

< 12.45 

< 12.46 

< 12.40 

< 12.40 

< 12.76 

< 12.76 

< 12.66 

< 12.66 

< 12.66 

< 12.26 

< 12.26 

< 12.10 

< 12.10 

< 12.10 

< 12.02 

< 12.02 

< 12.02 

< 12.02 

< 12.02 

< 12.02 

< 12.02 



15.06(0.02) 
14.90(0.05) 
15.09(0.03) 
14.87(0.07) 
14.98(0.05) 



13.31(0.02) 
13.33(0.02) 

< 12.44 

< 14.20 
14.75(0.05) 
14.76(0.10) 
14.29(0.18) 
14.58(0.05) 

< 14.20 
13.82(0.01) 
14.18(0.03) 
14.21(0.01) 
13.92(0.13) 
14.04(0.11) 
14.02(0.10) 
14.16(0.08) 
14.77(0.04) 
14.65(0.04) 
14.04(0.14) 
13.97(0.06) 
13.91(0.05) 
14.14(0.10) 
14.14(0.31) 
14.08(0.30) 

13.85(0.10)* 
14.05(0.38)* 
14.71(0.08) 
12.99(0.27) 
13.23(0.03) 
14.56(0.02) 
14.41(0.05) 
13.76(0.11) 
13.87(0.05) 
14.10(0.05) 
14.18(0.12) 
14.30(0.06) 
13.69(0.05) 
13.86(0.03) 
12.92(0.02) 
13.20(0.03) 
13.47(0.02) 
12.92(0.02) 
12.02(0.05) 
11.87(0.06) 
12.98(0.02) 
13.03(0.08) 
12.96(0.08) 



(0.02) 
(0.02) 
(0.09) 
(0.05) 
(0.02) 
(0.02) 
(0.02) 
(0.05) 
(0.28) 
(0.26) 
(0.03) 
(0.05) 
(0.05) 
(0.08) 
(0.08) 
(0.04) 
(0.03) 
(0.23) 
(0.03) 
(0.01) 
(0.03) 
(0.69) 
(0.02) 
(0.03) 
(0.02) 
(0.09) 
(0.06) 
(0.04) 
(0.16) 
(0.04) 
(0.09) 
(0.11) 
(0.06) 
(0.07) 
(0.18) 
(0.03) 
(0.10) 
(0.05) 
(0.05) 
^ 14.27 
14.86(0.02) 
sat/bid 
sat/bid 
sat/bid 
sat/bld 
sat/bld 
sat/bld 



15.101 
14.701 
12.29 
12.68 
13.33 
13.21 
13.61 
12.77 
12.501 
12.44 
14.78 
14.68 
14.26 
14.13 
14.22 
15.32 
15.09 
14.7^1 
15.39 
15.26 
14.89 
14.07 
14.90 
14.78 
15.22 
14.13 
12.93 
13.82 
13.90 
14.63 
13.23 
13.50 
13.31 
13.31 
12.72 
13.86 
12.99 
14.36 
14.61 



14.70(0.03) 


+ 0.03(0.06) 




10.0(0.3) 


Zn 


14.50(0.05) 


+ 0. 18(0.08) 


< 13.40 


5.5(0.4) 


Zn 


14.77(0.03) 


+ 0. 15(0.07) 


< 13.40 


8.4(0.5) 


Zn 


14.61(0.03) 


-0-69(0.05) 


12.82(0.04) 


8.7(0.7) 


Zn 


14.63(0.03) 


-0.13(0.04) 


12.79(0.04) 


8.2(0.2) 


Zn 




n 1 fi/fi nR"! 

— U. XDl^U.UO } 


< 12.92+ 


2.0(0.2) 




14.06(0.03) 


— 0.17(0.04) 


< 13.02"'" 


3.5(0.2) 


S 


< 14.14 


— 0.15(0.04) 


< 13.20 


3.5(0.2) 


Si 


< 14.14 


— 0.33(0.09) 


< 13.20 


6.9(0.4) 


Si 


< 14.14 


> +0.33 


< 13.20 


5.6(1. 1) 


Si 




> —0.61 


< 13.20 


7.2(0.8) 


Si 




— 0.48(0.06) 


< 13.00 


27.4(2.2) 


Si 




— 0.66(0.15) 


< 13.00 


14.9(3. 1) 


Si 




— 0.21(0.22) 


12.60 ± 0.08 


3.6(1.2) 


Si 




— 0.47(0.06) 


12.97 ± 0.06 


11.0(1.2) 


Si 




> —0.43 


< 12.60 


3.8(0.7) 


Si 




— 0.40(0.06) 


12.75±0.06 


7.2(1.0) 


Si 




— 0.66(0. 10) 


13.20ib0.04 


4.6(0.5) 


Si 




— 1.00(0.10) 


13.71±0.02 


11.8(1.1) 


Si 




— 0.18(0.13) 


12.56(0.24) 


4.7(0.3) 


S 




— 0.35(0.12) 


12. 14(0.72) 


6.8(0.8) 


S 




— 0.58(0.11) 


12.21(0.59) 


6.1(0.9) 


S 




— 0.37(0.09) 


12.70(0.17) 


3.7(0.3) 


S 




— 0.40(0.05) 


13.66(0.06) 


10.5(0.5) 


S 




— 0.38(0.05) 


13.26(0.07) 


6.3(0.3) 


S 




— 0.33(0.15) 


13.20(0.08) 


9.5(1.3) 


S 




— 1.13(0.09) 


12.61(0.20) 


2.5(0.6) 


S 




— 0.27(0.07) 


12.78(0.25) 


7.4(1. 1) 


S 




— 0.33(0.12) 


13.50(0.09) 


0.1(0. 1) 


S 




— 0.24(0.32) 


13.39(0.12) 


8.4(1.4) 


S 




— 0.24(0.32) 


13.38(0.12) 


5.5(1.0) 


S 




0.39(0.11) 


13.69(0.06) 


6.2(0.8) 


S 




0.08(0.38) 


13.26(0.10) 


3.9(0.4) 


S 




— 0.53(0.08) 


13.85(0.07) 


5.1(0.2) 


s 




— 0.84(0.82) 


11.83(0.29) 


0.3(0. 1) 


s 




— 0.19(0.04) 


12.75(0.05) 


9.8(0.9) 


s 




— 0.23(0.02) 


13. 12(0.03) 


5.3(0.2) 


s 




— 0.19(0.05) 


13.21(0.04) 


18.5(2. 1) 


s 




— 0.09(0.12) 


12.45(0.12) 


9.4(2.2) 


s 




— 0.31(0.06) 


12.38(0.07) 


3.8(0.6) 


s 




— 0.48(0.06) 


12.96(0.09) 


7.8(1.3) 


s 




— 0.26(0.12) 


13.16(0.12) 


6.2(1.4) 


s 




— 0.31(0.06) 


13.08(0.09) 


4.5(0.4) 


s 




— 0.19(0.07) 


12.75(0.03) 


5.1(0.5) 


s 




— 0.29(0.04) 


12.55(0.03) 


5.9(0.5) 


s 


< 13.37 


— 0.23(0.04) 


< 11.75 


7.5(0.2) 


Si 


< 13.37 


— 0.23(0.04) 


< 11.75 


5.4(0.2) 


Si 


< 13.37 


— 0.28(0.03) 


< 11.75 


16.0(0.6) 


Si 


< 13.67 


— 0.22(0.04) 


< 11.75 


4.4(0.2) 


Si 


< 13.67 


< +0.18 


< 11.75 


5.7(1.3) 


Si 


< 13.67 


< +0.33 


< 11.75 


2.0(1.3) 


Si 


< 13.24 


— 0.55(0.06) 


< 11.75 


10.9(0.5) 


Si 


< 13.24 


— 0.55(0.12) 


< 11.75 


6.9(0.9) 


Si 


< 13.24 


— 0.55(0.13) 


< 11.75 


3.6(0.3) 


Si 


14.26(0. 18) 


— 1.14(0.20) 


12.91(0.14) 


9.0(2.3) 


S 


14.43(0. 12) 


— 0.60(0.23) 


13.01(0.11) 


7.2(1.3) 


S 


13.47(0.27) 


— 0.51(0.34) 


< 12.70 


5.1(2.4) 


S 


14.42(0.05) 


— 0.26(0.07) 


< 12.70 


9.1(1.1) 


S 


13.35(0.27) 


— 0.40(0.30) 


< 12.70 


2.2(2.5) 


S 


14.65(0.01) 


— 0.20(0.03) 


< 13.50 


15.6(0.2) 


Zn 


14.23(0.01) 


— 0.05(0.04) 


< 13.08 


8.8(0.2) 


Zn 


< 13.05 


— 0.85(0.18) 


< 11.83 


4.5(0.4) 


Si 


< 13.05 


— 0.75(0.09) 


< 11.83 


10.3(0.3) 


Si 


13.19(0.26) 


— 0.49(0.02) 


< 11.83 


5.8(0.1) 


Si 


13.12(0.36) 


— 0.49(0.03) 


< 11.83 


9.4(0.2) 


Si 


13.32(0.19) 


— 0.45(0.02) 


< 11.83 


5.4(0.1) 


Si 


< 13.05 


— 0.23(0.05) 


b 


5.6(0.2) 


Si 


< 13.05 


— 0.15(0.28) 


b 


2.0(0.2) 


Si 


< 13.05 


— 0.42(0.26) 


b 


5.0(0.5) 


Si 


b 


— 0.27(0.06) 


< 13.66 


6.9(0.4) 


Zn 


b 


— 0.21(0.06) 


< 13.31 


6.5(0.6) 


Zn 


b 


> —0.47 


< 13.29 


6.2(0.3) 


Zn 




— 1.40(0.30) 


12.36(0.14) 


3.9(1.3) 


Si 




— 0.96(0.13) 


12.21(0.23) 


3.6(1.1) 


Si 




— 0.67(0.05) 


13,91(0.04) 


28.9(1,4) 


Si 




— 1. 16(0.06) 


13,74(0.04) 


6.0(0,4) 


Si 




— 1.10(0.31) 


13,00(1.37) 


6.1(1,5) 


Si 




— 1.21 (0.04) 


14,08(0. 18) 


5.6(0,0) 


Si 




— 1.07(0.02) 


13.87(0.02) 


7.2(0.0) 


Si 




— 0.92(0.05) 


13.96(0.03) 


8.1(0.5) 


Si 






13.14(0.06) 


10.3(0.0) 


Si 


14.56(0.03) 


— 0. 12(0.04) 


13.63(0.02) 


10.3(0.3) 


Zn 


14.54(0.04) 


+ 0.02(0.06) 


< 13.42 


14.0(0.8) 


Zn 


14.82(0.02) 


— . 05(0.03) 




14 7(0.3) 




< 13.85 




< 12.55 


4.3(0.7) 


Zn 


b 


-0.58(0.09) 


< 12.55 


5.0(1.3) 


Si 


b 


-0.40(0.04) 


< 12.55 


6.5(0.2) 


Si 




-0.22(0.17) 


b 


12.8(1.6) 


Si 




-0.34(0.04) 


b 


6.6(0.3) 


Si 


< 13.55 


-0.21(0.09) 


b 


6.7(0.8) 


Si 


< 13.55 


-0.26(0.11) 


b 


5.7(0.5) 


Si 




-0.29(0.14) 


< 13.00 


9.9(1.8) 


Si 




-0.27(0.13) 


< 13.00 


7.3(1.5) 


Si 




< +0.19 


< 13.00 


9.6(4.9) 


Si 


13.49(0.23) 


-0-36(0.04) 


< 12-53 


4.5(0.3) 


Si 


< 13.40 


-0-29(0-13) 


< 12-53 


5.8(6-1) 


Si 


13.97(0.05) 


-0-29(0-05) 


< 12-59 


5.7(0.3) 


s 


14.16(0.05) 


-0-25(0-05) 


< 12-59 


8.7(0.6) 


s 


< 13.56 


> -0-19 


< 12-59 


7.2(0.9) 


s 


14.50(0.02) 


-0-11(0-03) 


< 13-09 


8.5(0.2) 


Zn 


13.76(0.05) 


-0-50(0-05) 


< 12-52 


9.8(0.3) 


s 


13.49(0.07) 


-0-51(0-08) 


< 12-52 


5.1(0.5) 


s 


< 13.42 


> -0.44 


< 12.52 


3.2(0.4) 


s 


13.48(0.11) 


-0.09(0.11) 


< 12.52 


12.9(0.6) 


s 


13.85(0.04) 


-0.18(0.04) 


< 12.52 


4.9(0.2) 


S 


13.50(0.09) 


-0.55(0.10) 


< 12.52 


6.6(0.7) 


S 



The numbers in parentheses are la standard deviations. Strict upper limits are 5a- detection limits;" — ": no transition line from this ion has been observed 
"b" : blended; "sat/bld": most or all of the components in the observed transition lines from this ion are either saturated or blended, or both. 
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Figure 14. Voigt profile fits to systems with C II* absorption line and without H2 or C I absorption lines. In each panel, the histogram 
gives the observed spectrum and the continuous curve is the best fitted multicomponent Voigt profile. The vertical dotted lines mark the 
locations of individual components. 



— 2.0 ^ Zq ^ —1.5, are lower than what would be expected 
from CNM gas and are consistent with WNM (or low den- 
sity) gas. Interestingly these upper limits are lower than that 
seen in high latitude Galactic sightlines that are believed to 
be predominantly WNM gas. This means that the electron 
density (and therefore probably the total particle density) 
in these DLAs is probably quite small. 

6.3 Systems with H2 detection 

Systems with H2 detections (marked as filled circles in 
Fig. ESI have iV(C n*)/iV(C 11) consistent with CNM. We 
compute the allowed range of nn in these components us- 
ing A^(C 11*), A'^(S 11), rie from the C l excitation and 
T = Toi if available or T = 100 K (see EqlTJ. The resuhs 
are summarised in Column 6 of Table |1| This Table also 
gives upper limits on uh for components without H2 in sys- 
tems that show H2. It is to be remembered that we assume 
[C/S] in DLAs is [C/S]©. Realistically Carbon can be under- 
abundant by up to a factor of 2. In that case the density will 
be higher than what we quote in the table. From Table|3|and 
|l|it is clear that for the H2 components toward Q 0347—383, 
Q 0551-366 and Q 1232+082 the value of uh derived from 
both methods agree well. Such a comparison is not possible 
for the components toward Q 0013— 004 and Q 0405— 383 as 
C 11* is blended in the former case and C l is not detected 
in the latter case. In the case of Q 1444-1-014 the derived 
hydrogen density based on C 11* is lower than that derived 



Table 6. Systems without H2 detection 



QSO 


^aba 


log x(Al III) 


CNM^ 


WMN^ 


(om-3) 

Ionized^ 


Max* 


Q 0058 - 292 


2.671 




3.1 


1.2 


0.3 


<3 


Q 0112 - 306 


2.422 


-1.46 


24.4 


9.6 


2.3 


<15 


Q 0135 - 273 


2.799 


-1.86 


60.0 


23.4 


5.6 


<5 


Q 0405 - 443 


2.550 


-1.58 


7.3 


1.8 


0.4 


<4 


Q 0841 + 129 


2.374 


-1.16 


11.3 


2.8 


0.7 


<22 


Q 1157 + 014 


1.944 


-1.68 


16.3 


4.0 


1.0 


<3 


Q 1223 + 178 


2.465 


-1.35 


9.4 


2.3 


0.6 


<2 


1 T — 100 K ai 

3 T = 10^ K a 


,d „,/r 

,d Tlc/l 


iH = 0.001; 2 T 
^H = 0.1; * from 


= 8000 K a 
H2 cquilib 


r,d ,>e/n 


H = 0.01 





using C I fine-structure excitation. However, in this system, 
Ledoux et al. (2003) have found a 5 km/s shift between 
the C I absorption line and that of singly ionized species. 
In addition these components showing relative depletion of 
Si with respect to S, it is possible that we have over esti- 
mated A''(C 11). In summary, the nu estimates based on the 
two methods are approximately consistent with one another. 
The excitation of the fine-structure levels of C I and C 11 in 
the components with H2 detection are consistent with high 
density and low temperature CNM gas. 



6.4 Systems without H2 detection 

In this Section we focus our attention on the 7 systems in 
our sample that show C 11* without H2. These systems do 
not show detectable C I absorption lines except the high- 
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Figure 15. The average (over the whole profile) N{C n*)/N{C ll) 
ratio in DLAs (combining our sample to measures by Wolfe et al., 
2003a, 2003b) is plotted against log N(li l) (upper panel) and the 
metallicity Z (lower panel). The dotted lines show the expected 
range of values observed in the CNM gas with metallicity and 
dust-content ranging from that observed in the Galactic ISM to 
that of a typical DLA. The dashed lines give the corresponding 
range for the WNM gas. Filled circles and open squares are re- 
spectively for systems with and without H2 detections. In the 
top panel the filled triangles are the observations of Lehner et al. 
(2004) along high latitude Galactic sightlines with low H2 content. 
Most of the points with log N(H l) ^ 20.0 are from intermediate 
or high velocity clouds in the Galactic halo. 



metallicity system at Zabs = 2.1391 toward Tol 1037—270. 
Apart from the «abs= 1.943 system toward Q 1157-1-014 that 
show 21 cm absorption line (Wolfe et al. 1981) there is no 
independent constraint on T and nu- The identification of 
C II* at Zabs= 2.422 toward Q0112 + 029 and z^bs= 2.799 
toward Q 0135 — 273 is based on C ii*A1037 absorption line. 
These lines are well inside the Lyman-a forest and possi- 
ble contamination by intervening H l absorption cannot be 
ruled out. For the rest of the systems, the identification and 
estimation of the C 11* column density are secure. 

We detect Al lll absorption lines in the 6 (out of 7) 
systems for which our spectra cover the expected wavelength 
range of the redshifted Al ill transitions. We estimate the 
fraction of Al in Al ill using the observed metallicity and the 
observed A'^(A1 lll), assuming no depletion and solar relative 
abundances. Results are summarised in Table [5] From the 
Table, it can be seen that 1—7 % of Al is twice ionized. 
Using photoionization models from "CLOUDY" we derive 
a typical ionization parameter — 3 ^5 log U ^ —2 if the gas 
originates from a single slab irradiated by stellar spectrum 
with an effective black-body temperature of 30,000—40,000 
K (also see Fig. 1 in Izotov et al. 2001). This implies that 
the average We/nH ratio along the line of sight is typically 
in the range 0.3 to 0.9. Thus there are enough electrons in 
the cloud so that collisions with electrons are dominant in 
the excitation of C 11*. 

The average density, tih, is derived using Eq.|7|and as- 



suming three possible combinations of T and rie/nii- The 
results are summarised in Table |S| When we use no addi- 
tional constraints, the C 11* observations alone are consistent 
with the gas having high density and low temperature (see 
Column 4 of Table ISJ. 

The last column in the table gives the upper limit on 
TiH that will keep the equilibrium abundance of H2 below 
our detection limit (i.e A''(H2)^ 10^^ cm~^). This value is 
computed using simple formation equilibrium of optically 
thin H2 (Jura 1975) 



0.11/3(0)iV(H2) 



(8) 
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i?iV(HI) 

with, R and /3(0), respectively, the formation and photo- 
destruction rates of H2. In the case of the ISM, J? ~ 3x 10 
s-i cm"^ and (3(0) ~ 5 x 10"^° s"^ We use the ISM value 
of R scaled by the dust content measured in the systems. 
It can be seen from Table |S] that for a moderate radiation 
field (like the ISM mean field) tih derived using C 11* for 
the CNM like parameters is usually higher than the upper 
limit obtained based on the H2 equilibrium formation. In 
addition, the expected electron density for the T — 100 K 
gas (assuming a rio/riH ratio as seen in CNM) is higher than 
10~^ cm~"^. At such electron densities, C I should be de- 
tectable. This is inconsistent with the non-detection of C I 
in these systems. This problem of CNM gas producing very 
small amount of C I is already recognized in the literature 
(Liszt et al. 2002; Wolfe et al. 2003). 

We notice that the absence of C I and H2 in these sys- 
tems is consistent with the gas originating either from the 
WNM gas or from the ionized gas. As pointed out above, 
the strength of the Al III absorption lines seen in these sys- 
temsare consistent with the gas density being less than the 
one expected for the CNM. 

Thus, if one uses only C 11* absorption line then the 
results are consistent with these systems originating from 
CNM gas. However the absence of H2 and C I absorption 
lines together with the presence of Al ill following the pro- 
files of singly ionized gas is inconsistent with standard CNM 
solutions. Thus most of the DLAs without H2 are consistent 
with them originating from low density, high temperature 
and partially ionized gas. 



6.5 2abs= 1.944 toward Q 1157 + 014 

Some of our conclusions of the previous Section can be 
ascertained in the case of the 2:abs= 1.944 system toward 
Q 1157 -I- 014 as we have additional information on the ki- 
netic temperature based on 21 cm absorption (Wolfe et al. 
1981). The estimated spin temperature based on the recent 
measurement of N{R i) is T = 865±190 K (Kanekar & Chen- 
galur 2003). Fig. 1161 shows the velocity plot of selected ab- 
sorption lines in this system. The distribution of neutral, 
singly ionized and doubly ionized species in velocity space 
can be visualized easily. Using a Gaussian profile, we also 
show, in the bottom panel, the velocity range over which 
21 cm absorption is seen. Whereas the UV absorption lines 
spread over more than 100 km s~^, the 21 cm absorption 
originates from only a few of the components. There are 
components that do not possess 21cm absorption but show 
absorption due to UV transitions including that of C 11*. As 
expected, C 11* traces a wide range of physical conditions. 
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Figure 16. Velocity plot for the ^abs= 1-944 system toward 
Q 1157 + 014. The dotted Gaussian profile in the bottom panel is 
the reconstruction of the best fit Gaussian to the 21 cm observa- 
tions of Wolfe et al. (1981). This is just to illustrate the velocity 
range over which cold hydrogen gas is detected. 

Assuming the temperature of the gas that is producing the 
21 cm absorption feature to be 100 K (or 200 K), we can 
derive from the 21 cm observation the H I column density in 
the component lying along the line of sight: ~ 7 x 10^" cm~^ 
(or 1.4x10^^ cm"^). This is approximately 15% (30%) of the 
total A'^(H l) measured from the damped Lyman-a line. This 
means that 85% (70%) of A''(H l) along the line of sight is ei- 
ther warm or hot. From the Si 11AI8O8 profile, we notice that 
~60% of A''(Si 11) originate in the velocity space covered by 
the 21cm profile. We also notice that considerable fraction 
of Al III absorption originate from the velocity range covered 
by 21 cm absorption. Thus warm and ionized components 
seem to be co-spacial with the cold gas responsible for the 
21 cm absorption. In addition we also notice Al ill compo- 
nents with Si II and C 11* absorption well separated from the 
21 cm component. Thus part of the C 11* absorption seen 
here originate from the WNM or WIM. Therefore, the non- 
detection of H2 in this system can be easily explained as a 
consequence of most of the gas being at high temperature 
(and hence low density). The absence of H2 from the cold 
21 cm absorbing component could just be due to the low 
A'^(H l) associated with this component and the relatively 
low dust depletion. 



7 DISCUSSION AND CONCLUSIONS 

We have studied the physical conditions in damped Lyman- 
Q systems (DLAs) using a sample of 33 systems toward 26 
QSOs acquired for a recently completed survey of H2 in 
DLAs by Ledoux et al. (2003). We use standard techniques 
to estimate the physical conditions prevailing in the gas. 
In this Section, we discuss some of the results and related 
issues. 



7.1 High pressure of the H2 gas 

Our study shows that the H2 components in DLAs trace 
Cold gas (153±78 K) with relatively high pressure. The pres- 
sure in individual components (measured assuming a radia- 
tion field similar to our Galaxy) is in the range 824—30,000 
cm~"^ K, a large fraction of the components being at high 
pressure. 42%, 20%, and 8% of the components have pressure 
in excess of 3000 cm~^ K, 5000 cm"^ K and 10* cm'^* K, 
respectively. Based on the profiles of singly ionized species 
we note the H2 components arise in gas with a wide range 
of molecular content and ionization state much like what we 
see in the Galactic ISM. 

This is not unexpected. Indeed, in the framework of a 
galactic two-phase medium, the stable pressure range for 
the gas is 460=^ P/k (cm"^ K) ^ 1750 (Wolfe et al. 2003). 
Clearly the pressure we derive in the H2 components are 
much higher than this. From Table 3 of Wolfire et al. (1995) 
it can be seen that for a given metallicity an increase in the 
dust-to-gas ratio can lead to an increase in the allowed range 
of pressure, whereas an increase in the metallicity reduces 
the allowed range of pressure due to enhanced cooling. For 
conditions typical of DLAs, that is for metallicities of Z = 
— 1.0 and a dust-to-gas ratio ten times smaller than in the 
Galaxy, the stable pressure range is 1800-13000 cm~'^ s~^ 
(Wolfire et al. 1995). Note that in the absence of any confin- 
ing medium (or pressure equilibrium between different com- 
ponents) we expect such a high pressure gas to survive only 
for a short period of time (with a typical hydrodynamical 
time-scale of 10^ years). 

The pressure we infer depends very much on the in- 
tensity of the radiation field. A larger intensity implies and 
excess of UV pumping (on top of what we assume in our 
analysis) which, if taken into account, should reduce the hy- 
drogen density derived using the C I fine-structure lines. At 
the same time, the temperature of the gas will increase due 
to photo-heating. 

7.2 H2 content 

Ledoux et al. (2003) found that approximately 13—20% 
of DLAs show H2 absorption lines with most of the 
H2 components having column densities in the range 
16.0 ^ log iV(Il2)(cm"^) ^ 19.0. In the case of the Galactic 
ISM, only a minor fraction of the clouds fall in this range. 
This is expected because, above log A''(H2)— 16.0, self- 
shielding drastically decreases the photo-dissociation rate. 
On the contrary, in the LMC and SMC, a large fraction of 
lines of sight have 16.0 ^ log A''(H2) ^ 19.0 (see Tumlinson 
et al., 2002). Thus the trend noticed in DLAs could just be 
a generic feature of gas with low dust content and metallic- 
ity. It is however important to remember that the molecular 
fraction given in Ledoux et al. (2003) is an average over the 
whole line of sight. The actual molecular fraction in individ- 
ual components may be much larger. Thus the low values of 
7V(Il2) that are observed could just be a consequence of low 
A''(II l) in the corresponding individual components. Indeed, 
consistent models of DLAs (Srianand et al. 2005) require H I 
column densities much less than the total A'^(II l) measured 
in DLAs with H2. In addition, the absence of 21 cm absorp- 
tion can be reconciled if the II2 components have only part 
of the total H I (see below). 
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7.3 21 cm absorption 

We have shown that detecting H2 and C I absorption lines 
is an efficient way to trace the cold neutral gas in DLAs. 
H I 21 cm absorption line provides an independent way of 
detecting the CNM gas in DLAs. The detectability of 21 
cm absorption line depends only on the amount of cold gas 
along the line of sight and the covering factor of the ra- 
dio source. Thus, for a compact background source one can 
detect CNM gas with 21 cm without any bias from dust 
content or metallicity. At Zabs^ 2, seven systems have been 
searched for 21 cm absorption and none have been detected 
(Kanekar & Chengalur 2003). Assuming T = 200 K for the 
CNM gas, these authors estimated the filling factor of the 
CNM gas to be ^ 0.3. This is consistent with what we de- 
rive from our H2 survey. Over the redshift range covered by 
our survey there are three cases for which information on 21 
cm absorption and H2 content are available. To our surprise 
there seems to be no correlation between 21 cm absorption 
and H2 absorption. 

The 2abs= 2.811 system toward PKS 0528-255 show 
H2 absorption in two distinct components without any cor- 
responding 21 cm absorption (Carilli et al. 1996). The upper 
limit on r(21cm) gives iV(H !)< 5 x 10^° cm~^ if the kinetic 
temperature is similar to T(OPR) wo measure. Thus H2 and 
21 cm observations can be consistent with one another if no 
more than 20% of the total H I column density is associated 
with the H2 component. The 2abs= 1.944 system toward 
Q 1157-f014 (Wolfe, Briggs & Jauncy, 1981; and discussion 
above) and z^bs= 2.04 towards PKS 0458-020 (Briggs ot al., 
1989; Ge & Bechtold, 1999) show strong 21 cm absorption 
but no C I or H2 absorption. In these systems the absence of 
H2 and C I absorption could be either due to low density in 
the cold H I component or to the presence of higher ambient 
radiation field. 



7.4 C II* absorption 

As pointed out before, C 11* is detected in all the systems 
in which H2 is seen. In fact, C II* is also detected in the 
three systems that show 21 cm absorption discussed in the 
previous Section. C 11 being the dominant ion of Carbon in 
the neutral gas, it is natural to expect C 11* associated with 
both 21 cm and H2 absorption. However, C 11* is readily 
detected in warm neutral gas and even in ionized gas. The 
interpretation of the origin of the C 11* absorption is not 
as straightforward as in the case of H2 and C I. Thus, the 
nature of systems that do not show 21 cm absorption and/or 
H2 absorption is a matter of debate. The systems that show 
C II* in our sample are consistent with them originating from 
the CNM gas. However, the absence of C I and H2 (if we take 
the depletion as an indicator of the presence of dust) and the 
presence of Al ill are also consistent with C 11* absorption 
originating from the warm/partially ionized gas. Thus the 
frequency of occurrence of C II* provides a liberal upper 
limit on the CNM covering factor. A detailed investigation 
taking into account the constraints on the ionization state 
of the gas based on N 11, Fe ill or Al III will be important to 
derive the exact covering factor of CNM gas. 



7.5 Star-Formation Rate 

One of the main driver for the study of DLAs is to find out 
a way to recover the global star-formation history in a typ- 
ical, moderately star-forming environment. The importance 
of DLAs in the paradigm of hierarchical structure formation 
can be appreciated from the fact that the mass density of 
baryonic matter in DLAs at Zabs ~ 3 is similar to that of 
stars at present epochs (Wolfe, 1995). Studies of Lyman-a 
and UV continuum emission from galaxies associated with 
DLAs usually result in star formation rates (or upper limits) 
of a few Mq yr^i (Fynbo et al., 1999; Bunker et al, 1999; 
Kulkarni et al., 2001). 

Wolfe et al. (2003a, 2003b) have proposed a novel idea 
of using the C II* cooling rate to infer the SFR in DLA 
galaxies. The idea is that if one assumes thermal equilibrium 
then the cooling rate inferred from C II* should be equal to 
the heating rate (through UV photons. Cosmic rays etc.,) 
driven by the local star-formation activity. Their detailed 
study suggests the star formation rate density of DLAs at 
high- 2 could be as high as that inferred based on Lyman 
break Galaxies. We confirm the presence of C 11* in ~ 50% 
of the DLAs in our sample. However, it is clear from the 
above discussion that one needs to unveil the nature of the 
partially ionized gas in order to have a handle on the heating 
rates. 

In the local universe, star-formation is always related 
to molecular clouds. If DLAs are star-forming regions then 
the local star-formation has to be related to the mass of the 
molecular gas. Our survey shows that 13—20% of DLAs are 
associated with H2 in absorption. We have not detected CO 
in any of these systems and HD is detected in only one sys- 
tem (Varshalavich et al., 2000). Clearly the dark molecular 
clouds where stars form in our Galaxy are not seen along 
QSO lines of sight. The UV radiation field inferred from the 
H2 high-J excitation is similar to the Galactic mean field. 
Following Wolfe et al. (2003a,2003b) the SFR per unit co- 
moving volume for DLAs is, 

p, = Anco{z) < iiz) > 

■ , , , A .dN 

where, < ^(^;) > is the average SFR per unit area at redshift 
z and A, Ap and ^ are average physical cross-sectional, re- 
spectively, area, average projected area and number density 
of absorbers per unit absorption distance interval, /a is the 
fraction of DLAs in which the UV radiation field is similar 
to the Galactic UV background (i.e., 0.13—0.20). Here we 
use the fact that our H2 sample is a randomly chosen sub- 
sample of the whole population of DLAs and the presence 
of H2 is independent of A^(H l). 

For an Einstein-de Sitter cosmology, ^ = 3 x 10~^ for 
the mean redshift of our sample (Storrie-Lombardi & Wolfe, 
2000). Assuming Ho = 75 km s"iMpc"\ A/Ap = 2 and 
< ^{z) >= 4 X 10"^ M© yr"^ kpc"^ (typical for our Galaxy, 
see Kennicutt, 1998) we derive p. ^ 0.03 at 2abs= 2.5. This 
crude estimate already gives half the star formation rate 
density measured in Lyman break galax;ies (Steidel et al., 
1999). Recently, Hirashita & Ferrera (2005) have also arrived 
at similar conclusion. Thus it is of the utmost importance to 
understand the physics of the ISM in high-2 DLAs in order 
to derive the cosmic star-formation budget correctly. 
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